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SUMMARY 


Various  self-consistent  analyses  have  been  proposed  and  used  to 
approximately  evaluate  the  elastic  stiffness  and  elastic-plastic  behavior  of  metal 
matrix  composites.  Such  analyses  have  generally  relied  on  very  simple  theoretical 
models  for  the  matrix  inelastic  stress-strain  response.  This  was  perhaps  substantiated 
on  the  basis  of  a  lack  of  combined  stress  state  experiments.  Recent  biaxial  loading 
experiments  conducted  by  Dvorak  and  Bahei-El-Din  on  unidirectional,  fiber- 
reinforced  Boron-Aluminum  resulted  in  the  apparent  inapplicability  of  prior  self- 
consistent  approaches  based  on  (i)  the  measurement  of  "flats"  on  the  composite  yield 
surface  and  (ii)  nonnormality  of  the  plastic  strain  increment  with  respect  to  the  yield 
surface.  Dvorak  and  Bahei-El-Din  (1987)  introduced  a  new  theory,  the  so-called 
Bimodal  Theory  of  Fibrous  Composites,  in  response  to  these  apparent  deficiencies. 

Weng  (1988)  successfully  approximated  the  inelastic  behavior  of  spherical 
particle-reinforced  composites  utilizing  a  modified  self-consistent  model  called  the 
equivalent  inclusion-average  stress  (EIAS)  method.  Noting  the  overly  stiff  response 
of  the  basic  EIAS  model,  he  developed  the  "secant  modulus"  method  to  correct  for 
the  overconstraining  power  of  the  matrix. 

The  purpose  of  this  thesis  is  to  reexamine  the  problem  in  the  context  of  more 
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sophisticated  nonlinear  kinematic  hardening  rules  for  the  matrix.  An  EIAS  method 
which  incorporates  a  tangent  stiffness  formulation  based  on  incremental  plasticity  is 
proposed.  It  is  shown  that  this  method  is  comparable  to  Weng’s  secant  modulus 
method.  A  6  parameter  and  %  function  are  proposed  to  correct  for  the  constraining 
power  of  the  matrix  due  to  eigenstrain  accumulation  and  anisotropy  due  to  fiber 
reinforcement.  The  proposed  EIAS  method-tangent  stiffness  formulation  with  the 
6  parameter  and  %  function  produces  satisfactory  results  when  compared  to  existing 
experimental  data  for  the  Boron-Aluminum  system. 


CHAPTER  I 


INTRODUCTION 


Metal  matrix  composites  (MMCs)  represent  an  important  component  of 
today’s  materials  technology.  With  their  ductile  matrices  and  high-strength  fibers, 
MMCs  offer  high  strengthAveight  ratios,  superior  resistance  to  corrosive 
environments,  and  the  potential  for  high-temperature  applications.  Therefore,  in 
spite  of  their  high  production  costs  and  specialized  fabrication  techniques,  designers 
prefer  MMCs  over  the  more  widely  employed  resin  matrix  composites  (Pindera  and 
Lin,  1989)  for  many  applications. 

In  order  to  take  full  advantage  of  the  many  attractive  features  offered  by 
MMCs,  it  is  necessary  to  understand  and  be  able  to  accurately  model  their  response 
in  the  nonlinear  range.  While  researchers  have  enjoyed  reasonable  success  in 
modeling  the  inelastic  response  of  spherical  particle-reinforced  composites  (Tandon 
and  Weng,  1984;  Weng,  1988),  the  formulation  of  constitutive  relations  for  continuous 
fiber-reinforced  composites  has  been  frustrated  by  the  lack  of  reasonably  simple 
solutions  which  would  describe  the  fiber-matrix  interaction  and  its  effect  on  plastic 
deformation.  Progress  has  been  made  only  in  special  cases  where  certain  limitations 
are  placed  on  the  loading  conditions  or  constituent  materials  (Dvorak  and  Bahei-El- 


Din,  1982). 

Experimental  evidence  indicates  that  the  failure  strain  of  the  composite  will 
be  of  the  same  order  of  magnitude  as  that  of  the  elastic  fiber  (approximately  0.009 
for  boron,  for  example)  which  signifies  that  the  deformation  of  MMCs  is  limited  to 
the  range  of  small  strains.  In  this  situation,  plasticity  analysis  of  MMCs  is  best 
accomplished  using  the  micromechanics  approach,  which  not  only  produces  the 
effective  response  of  the  composite,  but  also  provides  an  indication  of  the  state  of 
stress  in  the  composite’s  constituent  phases.  Micromechanical  theories  calculate  the 
overall  response  of  the  composite  from  the  constituent  properties  and  fi-om  their 
mutual  constraints  which  are  given  by  microstructural  geometry  (Dvorak  and  Bahei- 
El-Din,  1982). 

Several  different  micromechanical  approaches  have  been  employed  previously  to 
model  the  inelastic  response  of  MMCs.  These  range  from  mixture  theories  for 
unidirectional  fiber-reinforced  composites  (c.f.  Murakami  and  Hegemier,  1986)  and 
self-consistent  schemes  (c.f.  Hill,  1965a;  1965b;  1965c)  to  finite  element  analyses  of 
a  repeating  unit  cell  (c.f.  Dvorak,  Roa,  and  Tarn,  1974).  The  most  recent 
developments  in  the  constitutive  theory  of  the  nonlinear  response  of  unidirectional 
MMCs,  based  on  the  micromechanics  approach,  include  the  models  proposed  by 
Dvorak  and  Bahei-El-Din  (1982,  1987)  and  Aboudi  (1986),  among  others  (Teply  and 
Dvorak,  1985  and  1988). 

Dvorak  and  Bahei-El-Din  have  made  significant  contributions  to  the  plasticity 
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analysis  of  fibrous  composites.  These  contributions  include  micromechanical  models 
which  derive  the  overall  response  from  uniform  local  fields  (averaging  models),  as 
well  as  models  which  approximate  the  actual  nonuniform  local  fields  and  arrive  at 
upper  and  lower  bound  solutions  (bounding  models)  (Bahei-El-Din  and  Dvorak, 
1989).  Perhaps  their  most  significant  contribution  in  this  area  came  as  a  result  of 
experiments  performed  on  unidirectionally  reinforced  B/Al  tubes  in  order  to 
characterize  the  elastic-plastic  behavior  of  composites  under  combined  stress  states. 
Based  on  their  observations,  they  developed  and  proposed  the  "bimodal  plasticity 
theory  of  fibrous  composites"  (Dvorak,  1987).  The  bimodal  theory  deals  primarily 
with  the  determination  of  composite  yield  surfaces  and  in  that  respect  may  be 
regarded  as  an  application  of  classical  irrcremental  plasticity  theory  to  heterogeneous 
composites. 

The  bimodal  theory  submits  that  the  elastic-plastic  response  of  MMCs 
reinforced  by  aligned,  continuous  fibers  can  be  described  in  terms  of  two  distinct 
modes.  In  the  matrix-dominated  mode,  the  composite  deforms  primarily  by  plastic 
slip  in  the  matrix,  on  planes  which  are  parallel  to  the  fiber  axis.  In  the  fiber- 
dominated  mode,  both  phases  deform  together  as  a  heterogeneous  medium  in  the 
elastic  and  plastic  range.  This  theory  admits  two  distinct  overall  deformation  modes 
which  may  exist  in  binary  elastic-plastic  fibrous  composite  systems  under  certain 
loading  conditions.  To  each  mode  there  corresponds  a  segment  of  the  composite 
yield  surface  which  reflects  the  onset  of  yielding  of  the  matrix  phase  in  that  particular 
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mode,  the  envelope  of  these  two  segments  comprises  the  composite  yield  surface. 
They  derived  constitutive  equations  for  the  matrix-dominated  mode  of  deformation 
while  the  response  of  the  fiber-dominated  mode  of  deformation  is  approximated  by 
a  self-consistent  model.  The  two  deformation  modes  give  different  branches  of  the 
composite  yield  surface  which  identify  the  state  of  stress  that  activates  a  particular 
mode,  and  indicate  the  conditions  for  mode  transition  in  a  given  composite  system. 
The  matrix-dominated  mode  is  found  to  exist  in  systems  reinforced  by  fibers  of  high 
longitudinal  shear  stiffness,  such  as  boron.  Systems  reinforced  with  more  compliant 
fibers,  such  as  graphite,  appear  to  deform  exclusively  in  the  fiber-dominated  mode 
(Dvorak  and  Bahei-El-Din,  1987). 

Dvorak  and  Bahei-El-Din  concede  that  some  observed  phenomena  cannot  be 
fully  explained  by  the  theory  and  appear  to  violate  some  of  the  classical  assumptions 
of  plasticity  theory  (eg.  Drucker’s  hypothesis).  For  instance,  they  observed  "flats"  in 
the  composite  yield  surfaces  and  non-normality  of  the  plastic  strain  rate  vector  to  the 
composite  yield  surface  (Dvorak,  1987;  Dvorak  and  Bahei-El-Din,  1988). 

Another  recent  work  which  stimulated  application  of  the  Equivalent  Inclusion- 
Average  Stress  (EIAS)  method  in  this  thesis  was  that  of  Weng  (1988).  Weng’s  work 
focussed  on  application  to  particle-reinforced  composites,  but  the  idea  of  relaxing 
constraint  in  the  EIAS  method  during  plastic  flow  is  especially  relevant  to  our  analysis 
of  fibrous-reinforced  composites. 

In  a  two-phase  composite  containing  spherical  inclusions,  two  distinct  modes 


4 


of  plastic  deformation  may  take  place.  The  first  one  is  that  the  matrix  will  remain 
elastic,  and  the  inclusions  are  the  plastically  deforming  phase.  The  second  one  is  that 
the  inclusions  will  remain  elastic  and  the  matrix  is  the  plastically  deforming  one.  As 
the  constraining  power 

imposed  on  the  deformation  of  the  inclusions  is  exercised  by  their  surrounding  matrix, 
these  two  types  of  composites  will  exhibit  different  constraining  effects.  In  the  former 
case,  as  the  matrix  remains  elastic,  its  constraint  power  will  remain  constant  and  be 
may  aptly  represented  by  its  elastic  stiffness  tensor.  Weng  terms  this  type  of 
composite,  with  plastically  deforming  inclusions,  a  composite  of  the  first  kind.  In  the 
second  case,  the  continuous  plastic  flow  of  the  matrix  phase  will  lead  to  a  weakening 
constraint  power  and,  therefore,  the  equation  describing  the  deformation  of  the 
inclusions  must  incorporate  such  an  effect.  This  type  of  composite,  with  the  matrix 
being  the  plastically  deforming  phase,  is  termed  a  composite  of  the  second  kind. 
Since  this  thesis  focuses  on  composites  of  the  second  kind,  we  will  direct  our 
attention  to  Weng’s  treatment  of  the  same. 

Weng’s  methods  for  analyzing  composite  plasticity  are  based  on  Mori-Tanaka’s 
concept  of  average  stress  in  the  matrix,  in  conjunction  with  Eshelby’s  solution  of  an 
ellipsoidal  inclusion.  Hence,  an  "interaction  strain"  is  introduced  in  addition  to 
Eshelby’s  transformation  strain  to  enhance  description  of  inclusion  interaction  effects 
at  higher  volume  fractions.  This  method  will  heretofore  be  referred  to  as  the 
equivalent  inclusion-average  stress  or  EIAS  method.  Weng  presents  two  solution 
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procedures  for  solving  the  problem  of  composite  plasticity  for  composites  of  the 
second  kind.  The  first  method  is  an  EIAS  method  with  elastic  constraint.  In  this 
method,  plastic  strain  in  the  matrix  is  treated  as  a  negative  eigenstrain  in  the 
comparison  material.  The  comparison  material  has  the  same  elastic  moduli  as  the 
matrix  material  and  remains  elastic  throughout  the  solution  procedure.  This  causes 
an  overconstraining  effect,  resulting  in  an  overly  stiff  composite  response  in  the 
plastic  regime  (Weng,  1988). 

Weng’s  second  solution  procedure  is  a  novel  EIAS  method  with  a  secant 
moduli  formulation.  Since  he  intended  to  apply  his  theory  for  monotonic, 
proportional  loading,  he  employed  the  deformation  theory  of  plasticity  instead  of  an 
incremental  theory  to  model  behavior  of  the  matrix.  With  a  ductile  matrix  material, 
there  is  decreasing  constraining  power  of  the  matrix  material  as  plastic  flow  increases. 
To  account  for  this  decreasing  constraint,  Weng  utilizes  the  "secant  modulus" 
approach  to  calculate  the  stiffness  of  the  comparison  material.  The  secant  modulus 
accounts  for  both  plasticity  and  the  weakening  constraint  power  of  the  matrix; 
therefore,  negative  eigenstrains  are  not  used  to  account  for  the  plastic  strain  in  the 
matrix  as  with  the  original  EIAS  method  with  elastic  constraint.  However,  c*, 
Eshelby’s  transformation  strain,  is  still  used  to  account  for  inclusion-matrix  misfit 
which  is  inherent  in  the  EIAS  approach.  Weng  notes  that  for  spherical  particle- 
reinforced  composites,  this  method  produces  much  better  results  than  the  EIAS 
method  with  elastic  constraint  (1988). 


The  purpose  of  this  thesis  is  to  develop  a  constitutive  model  which  predicts  the 
overall  response  of  the  composite  under  quite  generally  incrementally  applied 
uniform  macroscopic  stresses  or  strains  in  terms  of  the  constituent  properties  and  the 
geometry  of  the  microstructure.  Utilizing  Weng’s  analysis  of  the  composite  of  a 
second  kind  as  a  starting  point,  we  develop  three  material  models  for  analyzing  the 
plasticity  of  fibrous  composites.  Each  model  represents  a  refinement  of  the  previous 
in  order  to  achieve  a  better  approximation  of  experimental  data. 

First,  we  provide  a  detailed  development  of  the  EIAS  method  with  elastic 
constraint.  This  model  is  nearly  the  same  presented  by  Weng  (1988),  but  we  use  the 
Eshelby  tensor  of  elasticity  for  circular,  cylindrical  fibers  instead  of  the  tensor  for 
spherical  particle  reinforcement.  Also,  we  use  an  incremental  plasticity  theory  instead 
of  the  deformation  plasticity  theory  to  calculate  the  increment  of  plastic  matrix  strain. 
This  will  allow  us  to  subject  our  model  to  both  proportional  and  nonproportional 
loading  histories,  as  well  as  to  cyclic  loading  histories.  The  results  with  this  type  of 
model  as  applied  to  fibrous  composites  coincide  with  Weng’s  results  in  essence.  The 
model  performs  well  in  the  elastic  regime  but  is  overly  stiff  in  the  plastic  regime 
(Weng,  1988). 

Secondly,  we  develop  a  model  similar  to  Weng’s  EIAS  secant  modulus 
formulation;  however,  instead  of  using  the  secant  modulus  to  account  for  the 
weakening  constraint  power  of  the  matrix,  we  introduce  a  tangent  stiffness  method. 
This  method  uses  the  instantaneous  value  of  the  matrix  plastic  hardening  modulus  to 
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estimate  the  instantaneous  tangent  sti^ess  of  the  comparison  material  in  the  ElAS 
approach.  Our  use  of  an  incremental  plasticity  theory  to  compute  the  increment  of 
plastic  matrix  strain  instead  of  the  deformation  theory  of  plasticity  used  by  Weng 
motivates  the  introduction  of  the  tangent  stiffness  approach.  We  call  this  model  the 
BIAS  Method-Tangent  Stiffness  Formulation.  We  also  applied  this  method  to 
spherical  particle-reinforced  composites  and  compared  our  results  to  those  of  Weng’s 
using  his  secant  moduli  approach.  Our  results  are  comparable  to  Weng’s,  but  we  feel 
that  our  model  more  closely  approximates  experimental  data,  especially  at  lower  fiber 
volume  fractions. 

Thirdly,  in  light  of  the  unsatisfactory  results  achieved  with  the  above  two 
models  for  unidirectional  composites,  we  propose  a  theoretically-based  modification 
of  the  constraint  power  in  the  BIAS  formulation  which,  when  used  in  combination 
with  the  BIAS  method-tangent  stiffiiess  formulation,  produces  satisfactory  predictions 
for  the  inelastic  behavior  of  unidirectional,  continuous  fiber-reinforced  composites. 

Finally,  keeping  in  mind  the  experimental  observations  of  Dvorak  and  Bahei- 
Bl-Din  regarding  violation  of  normality  for  B/Al  continuous  fiber-reinforced 
composites,  we  offer  some  explanations,  based  on  our  theoretical  investigation. 
These  studies  lead  us  to  conclude  that  a  micromechanically-based  theory  is  essential 
for  describing  the  elastic-plastic  behavior  of  important  classes  of  fiber-reinforced 
composites.  Non-associativity  is  required  for  any  phenomenological  model,  in 
general,  owing  to  shear  localization  in  the  matrix  near  reinforcement  interfaces  in 
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addition  to  the  strong  tendency  for  the  plastic  strain  rate  in  the  matrix  material  to 
rotate  toward  the  direction  perpendicular  to  reinforcement  to  relax  constraint 
development. 

As  a  result  of  our  investigations,  we  propose  a  composite  material  model 
which  enables  fully  incremental  solutions  based  on  an  BIAS  method  with  eigenstrain 
relaxation  as  opposed  to  a  purely  rate  form  of  the  BIAS.  Our  original  contributions 
to  the  analysis  of  fibrous  MMCs  include  the  tangent  stiffness  formulation  for 
estimating  the  comparison  matrix  material  in  the  BIAS  method;  the  6  parameter 
which  relaxes  the  constraint  effect  of  the  comparison  material  by  reducing  the  rate 
of  eigenstrain  accumulation;  and  the  %  function  which  allows  the  inelastic  behavior 
of  the  composite  to  transition  to  matrix-  dominated  behavior  as  the  matrix  plastic 
strain  rate  vector  rotates  and  aligns  itself  with  the  reinforcement  direction. 
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CHAPTER  n 


EIAS  METHOD  WITH  ELASTIC  CONSTRAINT 

2.1  Eigenstrain  Terminology 

An  eigenstrain  is  a  collective  term  given  (Mura,  1987)  to  such  nonelastic 
strains  as  plastic  strains  and  misfit  strains.  It  is  a  "ficticious"  quantity  which  defines 
the  level  of  misfit  related  to  the  applied  loading  level,  mismatch  of  elastic  constants 
and  matrix  plastic  strain.  Eigenstress  is  a  generic  term  given  to  self-equilibrated, 
internal  stresses  caused  by  one  or  several  of  these  eigenstrains  in  bodies  which  are 
not  subject  to  any  externally  applied  body  forces  or  surface  tractions.  Eigenstress 
fields  develop  when  eigenstrains  are  incompatible. 

A  finite  subdomain  n,  with  a  prescribed  eigenstrain  Cy*,  in  a  homogeneous 
material  domain  D  with  no  prescribed  eigenstrain  is  called  an  inclusion  (see  Figure 
1).  The  elastic  moduli  of  n  and  D  are  assumed  to  be  the  same  when  inclusions  are 
under  consideration. 

If  n  has  different  elastic  moduli  than  D,  then  n  is  called  an  inhomogeneity. 
In  this  situation,  an  applied  stress  is  perturbed  by  the  presence  of  the  inhomogeneity, 
and  a  stress  perturbance  field  exists.  An  eigenstress  field  approximates  the  perturbed 
stress  field  by  considering  a  corresponding  fictitious  eigenstrain  Cy*  in  n  in  a 


homogeneous  material. 


2.2  Fundamental  Equations  of  Elasticity  with  Eigenstrains 

In  this  section,  the  equations  of  linear  elasticity  are  reviewed  with  particular 
reference  to  solving  eigenstrain  problems  (Mura,  1987;  Taya  and  Arsenault,  1989). 
These  problems  consist  of  solving  for  the  displacement  Uj,  strain  e^,  and  stress  Uy  at 
an  any  point  x(Xipc2PC3)  when  a  free  body  D  is  subjected  to  a  given  distribution  of 
eigenstrain  cy*.  A  body  is  free  when  there  are  no  externally  applied  surface  tractions 
or  body  forces. 


2.2.1  Hooke’s  Law 

For  infinitesimal  deformations,  the  total  strain  Cy  is  the  sum  of  the  elastic 
strain  ey  and  eigenstrain  Cy*, 


(M) 


Hooke’s  Law  relates  the  elastic  strain  to  the  stress  Uy  by 


Cyy  are  the  elastic  moduli.  In  D-n  where  the  Cy*  are  zero,  equation  (2-2)  becomes 


(2-3) 


In  this  work,  we  will  focus  on  isothermal  behavior  such  that  the  eigenstrains  are 
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2.2.2  Equilibrium  Conditions 

When  calculating  eigenstresses,  the  domain  D  must  be  assumed  free  of  any 
external  force  and  surface  constraints.  If  these  conditions  for  D  are  not  satisfied, 
superposition  of  the  eigenstress  field  of  D  on  appropriate  boundary  value  problem 
can  be  used  to  construct  the  applied  stress  field. 

The  equilibrium  equations  are 

<I„  =  0  (V  =  W3)  .  (M) 

The  boundary  conditions  for  free  external  surface  tractions  are 

=  0  ,  (2-5) 

where  nj  is  the  exterior  unit  normal  vector  on  the  boundary  of  the  domain  D.  In 
most  cases,  D  is  considered  an  infinitely  extended  body  and  the  boundary  conditions 
in  (2-5)  are  replaced  by  ajj(x)  0  as  x  -► «», 

23  Eshelby’s  Equivalent  Inclusion  Problem 

As  stated  earlier,  an  inclusion  is  defined  as  a  sub-domain  n  in  domain  D, 
where  eigenstrain  Cjj*(x)  is  given  in  n  and  is  zero  in  D-n.  The  elastic  moduli  in  n  and 
D-n  are  assumed  to  be  the  same.  The  remaining  domain  D-n  is  referred  to  as  the 
matrix. 

Eshelby  devised  a  method  to  solve  three-dimensional  elasticity  problems  for 
inclusions  (or  inhomogeneities)  of  ellipsoidal  shape  embedded  in  an  infinite  elastic 
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body  or  matrix.  Consider  first  an  inclusion  problem  where  an  ellipsoidal  domain 
(denoted  by  n)  is  subjected  to  uniform  non-elastic  strain  (or  eigenstrain)  (see 
Figure  1).  Such  an  inclusion  embedded  in  an  infinite  elastic  body  with  elastic  moduli 
Cyu  causes  stress  fields  within  and  outside  the  inclusion.  Eshelby  suggested  a  method 
for  obtaining  the  solution  of  the  stress  field  within  n  (Eshelby,  1957)  and  outside  n 
(Eshelby,  1959).  The  procedure  for  obtaining  the  stress  field  outside  n  is  very 
difficult,  but  that  for  the  stress  field  inside  n  is  simplified  because  it  is  uniform  for  a 
given  uniform  eigenstrain  Cy*.  In  the  absence  of  applied  loading,  the  perturbance 
stress  inside  Si  due  to  misfit  is  given  by  Eshelby  as 

«« = 

where  is  the  perturbance  strain  in  n  and  is  related  to  the  eigenstrain,  by 


'U  - 


Umn^  mn 


From  (2-6)  and  (2-7),  it  is  clear  that  the  perturbance  stress  inside  the  inclusion  can 
be  readily  calculated  for  a  given  eigenstrain,  Cy*. 

Syu  is  the  Eshelby  tensor  for  the  matrix  material  and  is  a  function  of  the 
geometry  of  the  ellipsoidal  inclusion  and  the  matrix  Poisson’s  ratio  if  the  matrix  is 
isotropic.  The  are  the  uniform  eigenstrains  in  the  inclusion  associated  with  the 
misfit.  Per  Eshelby,  uniformity  of  eigenstrains  requires  uniformity  of  the  remotely 
applied  stresses.  Sy^  is  given  in  terms  of  the  Green’s  function  G„,„pq  as 
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(2-8) 


1  2x 

-1  0 

where  the  integration  is  performed  over  the  unit  sphere. 

This  form  of  admits  any  form  of  anisotropy  of  the  matrix,  although  the 
Green’s  function  must  be  available  for  the  matrix.  The  above  equation  reduces  to 

The  G’s  are  presently  available  for  isotropic  and  transversely  isotropic  matrices.  For 
other  anisotropies,  they  may  be  determined  numerically. 

The  integration  for  Sp  is  performed  in  a  natural  coordinate  system  which 
defines  the  volume  of  an  ellipsoidal  inhomogeniety.  For  the  case  of  circular, 
cylindrical  fibers,  one  axis  of  the  ellipse  goes  to  infinity  while  the  other  two  are  equal. 
For  a  more  detailed  development  of  the  Eshelby  tensor,  see  Appendix  A. 

A  practical  and  important  application  of  Eshelby’s  equivalent  inclusion  method 
is  the  case  of  an  "inhomogeneity  problem."  When  the  elastic  moduli  of  an  ellipsoidal 
inclusion  differ  fi'om  those  of  the  surrounding  matrix,  the  inclusion  is  called  an 
ellipsoidal  inhomogeneity.  Fibers  and  precipitates  are  examples  of  inhomogeneities. 
An  elastic  material  containing  inhomogeneities  is  free  from  any  stress  field  unless  a 
load  is  applied.  On  the  other  hand,  a  material  containing  inclusions  is  subjected  to 
a  self-equilibrated,  internal  stress  (eigenstress)  field,  even  if  it  is  free  from  all  external 
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load  is  applied.  On  the  other  hand,  a  material  containing  inclusions  is  subjected  to 
a  self-equilibrated,  internal  stress  (eigenstress)  field,  even  if  it  is  free  from  all  external 
tractions.  Eshelby  first  pointed  out  that  the  stress  perturbance  in  an  inclusion  relative 
to  the  applied  stress  due  to  the  presence  of  an  inhomogeneity  can  be  simulated  by 
an  eigenstress  caused  by  an  inclusion  when  the  eigenstrain  is  chosen  properly.  This 
eigenstress  is  called  the  perturbance  stress,  and  the  equivalency  will  be  called  the 
equivalent  inclusion  method. 

When  a  composite  consisting  of  a  matrix  with  stiffness  tensor  Qjkl  and  a  single 
inhomogeneity  with  a  stiffness  tensor  Qjkl  is  subjected  to  the  uniform  applied  stress 
a■^°,  the  actual  (or  total)  stress  field  in  the  inhomogeneity  (see  Figure  2(a))  is 
given  by 

where  and  eyP'  are  the  stress  and  strain  perturbances  introduced  by  the  existence 
of  the  inhomogeniety,  and  is  the  strain  of  the  matrix  (without  the  inhomogeneity) 
corresponding  to 

This  inhomogeneity  problem  can  be  reduced  to  the  inclusion  problem  insofar 
as  the  disturbances  of  the  stress  and  strain  are  concerned.  In  the  equivalent  inclusion 
method,  we  must  introduce  an  eigenstrain  in  n  in  order  to  reproduce  the  effect  of  the 
presence  of  the  inclusion.  The  composite  is  treated  as  though  it  consists  of  pure 
matrix  material  (elastic  constants  Qjij).  There  are  no  eigenstrains  introduced  in  D-n. 


15 


Hence,  in  n,  we  have 


oJ+oy=Cja(ea+e^  actual 

^Quivalent  inclusion  problem 


(2-11) 


where  Cy*  is  the  eigenstrain.  Equation  (2-11)  is  the  governing  equation  of  misfit  (see 
Figure  2(b))  for  a  single  inhomogeneity. 

It  is  obvious  from  equations  (2-10)  and  (2-11)  that  the  perturbance  stress  in 
the  inclusion  is  given  by  equation  (2-6).  Then,  the  perturbance  strain  is  related 
to  the  unknown  eigenstrain  Cy*  by  equation  (2-7).  The  perturbance  stress  in  n  is 


“  ®«)- 


(2-12) 


Note  that  ayP'  -♦  0  and  Cy^'  -*  0  as  jxjj  -♦  »,  such  that  is  self-equilibrated,  i.e. 


V'=o. 


If  the  matrix  has  uniform  plastic  strain,  CyP®,  then  we  may  treat  CyP®  as  a 
negative  eigenstrain,  i.e.. 


=  -^ii  ^  «v- 


(2-13) 


In  the  context  of  the  elastic  Eshelby  solution,  the  matrix  plastic  strain  must  be 
enforced  homogeneously  in  D  and  then  treated  as  a  negative  eigenstrain  in  n  (Weng, 
1988;  Mura,  1987).  In  metal-matrix  composites  for  which  the  matrix  is  elastic-plastic 
and  the  fibers  are  elastic,  plastic  strain  occurs  in  the  matrix.  We  must  then  modify 
(2-11)  as  follows  (Weng,  1988): 
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,»•  _  p**\ 

' «  ®  « /» 


where  it  is  understood  that  CuP“  is  the  matrix  plastic  strain.  The  left  side  of  (2-14) 
represents  the  actual  stress  in  the  inclusion  assuming  the  inclusion  deforms 
compatibly  with  the  matrix  plastic  strain  (i.e.  Ui'*^=Uj  ).  The  right  side  is  the  associated 
equivalent  inclusion  problem.  The  perturbance  strain  in  this  case  is 


(2-15) 


It  is  noted  here  that  the  equivalent  inclusion  problem  employs  the  assumption 
of  elastic  constraint  of  the  matrix  on  the  fiber  which  means  that  the  problem  treats 
the  matrix  as  if  no  plastic  strain  had  occurred.  This  may  be  reasonable  in  the  sense 
that  after  plastic  straining,  misfit  accommodation  is  by  virtue  of  elastic  unloading 
processes,  i.e.  within  the  yield  surface,  at  each  successive  increment  of  deformation. 
Weng  (1988)  pointed  out  that  this  procedure  overestimates  the  constraint  power  of 
the  matrix  at  significant  levels  of  matrix  plasticity. 


2.4  Interaction  of  Inclusions 

The  development  of  the  inhomogeneous  inclusion  problem  was  based  on 
Eshelby’s  equivalent  inclusion  concept  for  a  single  inclusion  in  an  infinitely  extended 
matrix  material.  This  approach  needs  to  be  modified  when  it  is  used  to  analyze 
fibrous  composite  systems.  In  a  fibrous  composite,  each  fiber  is  an  inclusion. 


Depending  on  the  volume  fraction,  the  fibers,  or  inclusions,  tend  to  interact  causing 
additional  stresses  in  n  which  augment  the  perturbance  stress,  These  additional 
stresses  are  called  interaction  stresses,  in  general,  and  are  represented  by  a^^. 

Strictly  speaking,  we  cannot  just  assume  for  finite  D  that  on  the 

boundary  of  D  as  is  the  case  with  a  single  inclusion  in  an  infinite  matrix.  In  general, 
we  may  think  of  Uy*  as  an  image  stress,  analogous  to  image  stresses  due  to  finite 
boundaries  in  dislocation  mechanics,  for  example.  Then,  equation  (2-14)  becomes, 
in  n. 


Cjt^ea+e^+c^+ei)  = 

(Z-16) 

Here,  Cy'  and  are  regarded  as  an  averaged  image  or  interaction  strain  and  stress, 
respectively,  over  the  complete  volume.  Then  we  must  specify  how  Cy’  is  to  be 
calculated. 


2.4.1  The  Mori-Tanaka  Method 

The  Mori-Tanaka  Method  assumes  that  the  inclusions  are  randomly 
distributed.  Then,  per  a  rule  of  mixtures  type  formulation,  in  the  absence  of  remote 
applied  loading, 

+  (1  -  =  0- 

Here,  (uy)^  =  the  average  internal  (self-equilibrated)  stress  perturbation  in  the 
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inclusion; 

(aij>nj  s  the  average  internal  (self-equilibrated)  stress  perturbation  in  the  matrix; 
and  S  =  the  volume  fraction  of  inclusions. 

The  total  stress  inside  an  inclusion  is 

=  oj  +  (o..)^  (2-18) 

where 

<  -  Cw/'S  -  *  eD  (2-l») 

and 

eS  -  -  O-  <«*) 

Here,  ajj*"  is  the  stress  corresponding  to  the  infinite  body  solution,  i.e.  a  single 
inclusion  in  a  matrix  of  infinite  extent. 

Since  the  inclusion  can  be  placed  anywhere  in  the  matrix,  the  average  internal 
stress  in  an  inclusion  is 

<Otf>/  =  o? 

Substituting  for  (ay)^  in  (2-17)  and  simplifying,  we  have 

*  (1  -  =  0 

f<  -  (2-22) 
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Utilizing  this  result  in  combination  with 


<Oy>/  =  (2-23) 

which  was  obtained  by  dividing  through  by  /  in  (2-17),  we  have 

<(.„>,  -  (1  -  /)of.  (2-M) 

Let  us  define  the  average  interaction  stress,  i.e. 

o".  -  (o,>  .  (2-25) 

!f  '  (f' m 

Therefore,  according  to  the  Mori-Tanaka  Method,  the  interaction  (or  image)  stress 
is 

aj.  =  -/oj  .  (2-26) 


2.4.2  The  Image  Stress  Approach 

Mura  (1987)  presents  the  concept  of  an  image  stress,  where  the  internal 
stress  is  equal  to  the  sum  of  the  image  stress  and  the  perturbance  stress,  ayP',  in  n, 
i.e. 


/  pt 

+  oj  . 


(2-27) 


Consider  a  finite  material  domain  D  in  which  are  embedded  n  number  of 
ellipsoidal  inclusions.  Also,  assume  there  exists  some  finite  ellipsoidal  domain,  D*, 
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which  encloses  the  inclusion,  n,  and  has  dimensions  on  the  order  of  the  fiber  spacing 
(see  Figure  3(a)).  The  domain,  D*,  is  of  the  same  ellipsoidal  character  as  n  in  terms 
of  orientation  and  shape  (see  Figure  3(b)).  It  can  be  shown  that  aij^*  by  itself  (i.e. 
a  single  inclusion  in  an  infinite  body)  does  not  satisfy  the  traction  boundary  condition 
aij‘’*nj=0  on  tfD*.  Hence,  the  image  stress  may  be  viewed  as  a  corrective  term  to 
satisfy  equilibrium  in  D*  and  the  traction  boundary  conditions  on  5D*,  i.e. 


Ogj  =  0  in  D\ 

(oj  +  af)nj  =  0  on  6D\ 


(2-28) 


Recall  that  ajjjP‘=0  from  the  equivalent  inclusion  problem.  This  ensures  some 
measure  of  consideration  of  the  finite  distance  between  fibers. 

The  average  value  of  over  the  volume  D*  can  be  shown  to  be 


<<’'»>  =  --tL  <*•») 

" D*  '  p*  “ 

by  Tanaka-Mori’s  Theorem  (Mura,  1987).  V^.  is  the  volume  of  the  domain  D*,  and 
the  integration  is  performed  over  the  volume  of  the  inclusion  n. 

For  uniform  eigenstrains  in  n,  is  constant  in  n  and  hence 


-ro 


pt 

a 


(2-30) 


I  ^  i 

Here,  n  is  the  volume  of  the  inhomegeneity,  and  n/Vjj.  =  /*,  the  fiber  volume 
I 
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fraction  in  D*.  In  our  analysis,  we  then  assume  that 


(oj)  =  aj  .  (2-31) 

in  both  phases.  Likewise,  for  strain, 

<'«>  = 

and 

<ej>  -  -/-K  -  '■*  "  *r). 

where  it  is  understood  that  the  brackets  refer  to  the  average  over  the  volume  D*. 

Utilizing  this  result  and  attributing  to  the  bracketed  quantities  which  follow  the 
definition  of  average  over  the  entire  volume  of  the  finite  domain  D,  we  are  able  to 
determine  an  expression  for  the  image  stress  in  terms  of  the  perturbance  stress  for 
the  entire  composite.  The  image  strain,  is  related  to  the  image  stress,  a  ^,  by 

ei  >  .  (2-34) 

From  (2-19),  (2-33),  and  we  have 

oj  •  -/of  ■  c4-/(.S  -  eL  *  e{D]  .  (2-35) 

The  quantity  in  brackets  on  the  right  side  of  the  equation  is  as  can  be  seen  from 
(2-33).  Note  that  we  are  now  using  /,  the  composite  fiber  volume  fraction,  instead 
of  f.  Recognizing  this  and  solving  for  Cu‘,  we  have 


(2-36) 


We  can  clearly  see  that  (2-34)  and  (2-36)  are  both  expressions  for  Equating  the 
right  sides  of  these  equations  and  multiplying  both  sides  by  Gyu,  we  see  that 

cj  .  -/oj  (2-37) 

and  therefore  the  average  interaction  strain  in  (2-33)  is  consistent  with  (2-30).  One 
notes  that  the  above  expression  for  the  image  stress  is  identical  to  that  obtained  by 
the  Mori-Tanaka  method,  which  may  be  expected  for  a  two-phase  composite  with 
uniform  eigenstrains  and  ellipsoidal  inclusions  as  pointed  out  by  Norris  (1989). 


2.5  Average  Stress  and  Strain  in  the  Composite 

Stress-strain  analysis  usually  requires  some  assumption  about  the  relationship 
of  the  applied  stress  to  the  actual  stress  in  the  composite.  A  usual  approach  is  to 
take  the  average  of  all  the  stresses  over  the  entire  volume  of  the  composite,  i.e. 


=  —f  ojiV 

^  vJyD  ^ 


(2-38) 


Since  the  various  stresses  are  assumed  uniform  in  the  BIAS  approach,  we  have 


=  (1  -  /)(oJ  +  oj)  +  /(oj  1-  oj  +  of).  (2-39) 

Hence,  the  average  stress  in  the  composite  is  equal  to  the  matrix  volume  fraction 
times  the  average  matrix  stress  plus  the  fiber  volume  fraction  times  the  average  fiber 
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stress.  Simplifying  (2-39),  we  have 


—  O  I  X  pt 

=  0(,  +  °ij  f°ij 


In  section  2.4,  we  showed  that 


Therefore, 


(^40) 


(2^1) 


Oy  =  Oy  • 


(242) 


Hence,  the  applied  stress,  ay°,  is  equal  to  the  average  stress.  This  should  be 
intuitively  obvious  based  on  the  foundations  of  the  EIAS  approach. 

Applying  a  similar  volume  averaging  approach  to  determine  the  average  strain 
in  the  composite,  we  have 


8^  =  ®r  (1-^(8 J+cJ)  +  /(8j+8j+E^).  (243) 

Therefore,  the  average  strain  in  the  composite  is  equal  to  the  sum  of  the  matrix 
plastic  strain  (corresponding  to  a  homogeneous  stress-free  deformation),  the  matrix 
volume  fraction  times  the  average  matrix  strain,  and  the  fiber  volume  fraction  times 
the  average  fiber  strain.  After  simplification,  we  have 


(244) 


—  pm  o  I  X  Pt 

®</  =  ®v  '®^ 

=  (1  -  /)8j"  +  cj  +  /e*  . 

Note  that  we  have  foregone  any  consideration  of  variational  principles  for 
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defining  the  strain  of  the  aggregate.  These  considerations  are  not  essential  when 
dealing  with  a  work-hardening  matrix. 

2.6  Final  Form  of  the  ElAS  Method  with  Elastic  Constraint 

At  this  point  we  are  ready  to  combine  the  concepts  of  the  equivalent  inclusion 
problem  with  those  of  image  and  average  stresses  to  produce  the  final  constitutive 
relationships  for  our  model  based  on  the  treatment  of  matrix  plastic  strain  as  a 
negative  eigenstrain  in  the  fiber. 

For  better  description  of  finite  fiber  volume  interaction  effects,  the  image 
stress  is  included  on  both  sides  of  (2-14).  Repeating  (2-16),  in  n, 

^  (245) 

=  +  eS  t-  +  cifc,  -  e*«)  . 

The  left  side  of  (245)  represents  the  actual  stress  in  the  fiber  whereas  the 
right  side  represents  the  corresponding  equivalent  inclusion  problem  for  the  same 
stress. 

2.7  Rate  Form  of  the  EIAS  Approach 

We  may  consider  the  solution  of  the  elastoplastic  problem  as  the  sum  of  the 
solutions  of  a  sequence  of  elastically  constrained  plastic  deformations,  subjecting  each 
increment  to  the  equivalent  inclusion  solution. 
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Adopting  a  rate  form  of  equation  (2-45),  we  have 


♦  tS  *  ej-  •*  - 


where 


«  O  •  •  Bin  •  I  *  ^  V 


®«  “  ^kbnK^^mn  » 

««  =  -/(c«  -  ®W  +  CtT)  » 


and  iu  =  CjJdJ  . 


(246) 


(247) 


The  matrix  plastic  strain  rate,  c  is  given  by  the  incremental  elastoplastic  flow  rule 


cf  =  0  if  F  <0  or  F  ^  0  and  —a^  <  0, 

do; 

ej"  =  if  F  =  0  and  ^  0- 


(248) 


Here,  F  s  the  matrix  yield  function,  and  =  0^°  +  (7^*  is  the  average  matrix  stress 
rate.  The  unit  normal  vector  to  the  yield  surface  is  given  by 


-  JLlJLJLjz 


daJdOjudOt, 


(249) 


The  quantity  h  is  the  plastic  hardening  modulus  of  the  matrix.  Hence,  we  may  restate 
(2-48)  as 


(2-50) 


if  ^  0  and  F  =  0, 

=  0  otherwise. 

The  quantity  enclosed  in  brackets  is  the  matrix  stress  rate,  <7^™.  Appendix  B  outlines 
the  two  incremental  plasticity  theories  used  to  calculate  the  matrix  plastic  strain 
increment  in  this  thesis. 

Substituting  (2-47)  into  (2-46),  we  have 

-  O  - 

-  -  o  - 1;  *  *rl}  - 

(2-51) 

cJcLoL  *  -  O  *  'S' 

“  'SI)  ■  'U  +  'S']  ■  '«}  • 

In  this  equation,  a^°  is  the  applied  stress  rate  (or  increment)  and  the  unknowns  are 
eu*  and  c/*". 

2.8  Solution  Procedure  for  the  ElAS  Method  with  Elastic  Constraint 

Equation  (2-51)  above  is  the  governing  equation  for  solution  of  the  EIAS 
problem  with  full  elastic  constraint.  Since  an  incremental  plasticity  theory  is  used  to 
calculate  the  increment  of  matrix  plastic  strain,  an  incremental  form  of  (2-51)  is 
adopted  for  solution  of  the  problem.  Equation  (2-51)  is  rewritten  for  ease  of 
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programming  as  follows: 


where  the  coefficient  matrices  and  Djj^„  are 


(2-52) 


^i^mn  ^  ~^I^Umn  *  *  ^(fk^Uiat  ~  ^^tfl^Umn 

*  ~  ^ijmn  ^iimn  ~ 

=  (1  -  fX^^iimn  ~  ^^ak^Umnj 
^aum  ~  ^vk^Umii  ~  ^^yl^klmn  ~  ^Vk^Umn  *  f^yU^kk 
f^ijnui  '*'  "  ^(fmn 

=  (1  -  fx-^c^^  -  c^)  -  /c;^ . 


where  ACjj„,n  —  Cymn  -  CjjjQ„. 


(2-53) 


These  matrices  are  calculated  based  on  the  material  properties  of  the 
particular  composite  system  of  interest,  and  they  only  have  to  be  calculated  once  for 
a  given  loading  history.  Therefore,  the  subroutine  which  performs  these  calculations 
(called  COEFMAT  for  material  coefficients)  is  located  outside  the  calculation  "loop", 
which  significantly  reduces  calculation  time. 

The  program  was  written  for  "load  controlled"  conditions  as  the  applied  stress 
increment,  Aay®,  is  the  input  to  the  program.  A  separate  program  capable  of 
producing  proportional  as  well  as  nonproportional  loading  histories  was  written  to 


create  the  loading  history  input  file. 

As  stated  above,  the  two  unknowns  in  the  problem  are  the  increment  of 
eigenstrain,  ac^*,  and  the  increment  of  matrix  plastic  strain,  ACyP™.  The  plastic  strain 
increment,  ACyP™  is  calculated  based  on  the  current  stress  state  (and  increment)  using 
an  incremental  plasticity  theory  outlined  in  Appendix  B.  With  Aay®  and 
known,  ac^*  can  be  easily  calculated  by  multiplying  the  left  side  of  (2-52)  by  the 
inverse  of  the  coefficient  matrix  Djjn,„  (i.e.,  Djj^n'^).  Then,  the  complete  state  of  stress 
and  strain  for  the  composite  and  its  constituents  can  be  calculated  for  that  particular 
instant  in  the  loading  history.  Updating  of  the  stress  and  strain  was  accomplished  in 
the  program  by  a  subroutine  called  UPDATE.  Once  the  update  was  performed  and 
values  were  written  to  a  data  file,  the  program  reads  the  next  applied  stress 
increment  from  the  input  file  and  the  cycle  continues  until  the  loading  history  is 
exhausted. 

This  program  is  particularly  efficient  in  that  no  iteration  is  required.  This  is 
accomplished  as  follows.  From  (2-50),  we  have  the  matrix  flow  rule 


= 


Taking  the  scalar  product  of  both  sides  of  (2-54)  with  n^j,  we  have 

i  +  eg]n„. 

From  the  EIAS  method,  equation  (2-52), 


(2-54) 


(2-55) 
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(2-56) 


f{^Umn^mtip9  ^Up^r^rsafi  tu^U _  (2-59) 

This  noninterative  form  makes  the  program  especially  practical  for  application  in 
finite  element  type  programs,  for  example.  Readers  are  referred  to  Appendix  C  for 
a  flowchart  of  this  program,  called  BALNIT  (for  Boron/ALuminum-No  ITeration). 


CHAPTER  m 


EIAS  METHOD  -  TANGENT  STIFFNESS  FORMULATION 


3.1  EIAS  Method  with  a  Tangent  Stiffiaess  Formulation 

Based  on  his  work  with  spherical  inclusions  in  a  plastically  deforming  matrix, 
Weng  (1988)  concluded  that  the  constraint  power  of  the  EIAS  method  with  elastic 
constraint  renders  it  much  too  stiff  in  the  elastic-plastic  range.  Weng  introduced  a 
modified  Eshelby  formulation  based  on  the  secant  modulus  for  the  deformation 
theory  of  plasticity  in  order  to  weaken  the  constraint  power  of  the  matrix  in  the  EIAS 
method.  For  incremental  elasto-plasticity  and  general  loading  paths,  we  caimot  use 
the  secant  modulus  approach  directly.  However,  it  is  possible  to  construct  an 
analogous  approximate  Eshelby  relation  based  on  the  tangent  modulus.  It  is 
important  to  emphasize  that  the  EIAS  method  (Weng,  1988;  Mura,  1987)  based  on 
elastic  constraint,  which  treats  plastic  strain  as  an  eigenstrain,  is  not  generally 
applicable  to  significant  matrix  inelastic  deformation  of  reinforced  materials.  Even 
though  the  method  may  be  applied  unambiguously  to  the  case  of  viscoplastic  matrix 
behavior,  it  is  nonetheless  overly  stiff. 


To  introduce  a  tangent  stiffness  formulation,  we  note  that  the  matrix  plastic 
strain  rate  is  not  explicitly  required  as  in  the  previous  case.  The  plastic  hardening 


modulus  h,  based  on  the  plastic  internal  state  of  the  matrix  material,  is  uniquely 
related  to  the  tangent  stiffness, 

For  deviatoric  matrix  plasticity,  we  may  write 


Solving  for  the  matrix  strain  rate  during  plastic  flow. 


(3-1) 


= 


t  T  \“^  •  •• 


^kl 


(3-2) 


We  may  write  the  modified  Eshelby  rate-type  equation  based  on  tangent  stiffiiess  as: 

*  ®«)  ®  ~  ®tf)- 

Since  the  effect  of  matrix  plastic  deformation  is  introduced  through  the  tangent 
stiffness,  we  note  that 


®  iU  ~  mn* 


.  / 


Cy  =  -/(®?  -  «i). 

and  =  (clSf^Ou 


(3-4) 


where  the  Mori-Tanaka  theorem  (Mura,  1987)  is  employed  to  arrive  at  the  form  of 
the  interaction  strain  rate  Cy*.  It  is  important  to  note  that  a  tangent  stiffness  method 
must  be  employed  in  the  reference  material  relation  between  Cy®  and  ay°;  strictly 
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speaking,  should  be  based  on  the  elastic-plastic  state  of  a  comparison  matrix 
material  at  stress  0^°,  but  since  and  are  no  longer  uniquely  related  (as  is  the 
case  for  deformation  theory),  a  formulation  for  will  be  introduced  later. 

The  Eshelby  tensor  in  this  case  is  based  on  the  effective  Poisson’s  ratio, 
V,  of  the  matrix  which  may  be  assumed  of  the  form 


V  =  V  + 


(0.5  -  v) 


(3-5) 


where  v  is  the  elastic  Poisson’s  ratio  of  the  matrix.  Hence,  v  varies  continuously 
from  V  for  h  -»  00  to  0.5  at  h  =  0.  This  form  is  based  on  the  dependence  of  the 
tangent  stiffness  on  the  plastic  modulus  h.  The  tangent  stiffness  formulation  is  rather 
insensitive  to  the  choice  of  the  form  of  V  so  long  as  V  is  allowed  to  vary  continuously 
between  fully  elastic  and  fully  plastic  limits. 

The  average  stress  rate  in  the  matrix  is  given  by 


(3-6) 


where 

oj  =  -/dj  =  -/C^ycS  -  cy  (3-7) 

in  the  tangent  stiffness  formulation.  The  average  composite  stress  and  strain  rate  are 
given  by 
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Og  =  (1  -  f)iol  +  a‘^  +  /(dj  +  dj  +  djf) 

=  a^j  *  Off  +  fa^  = 

‘e^  =  (1  -  ml  +  ej)  +  /(ej  +  ej  +  ej) 


=  Sy  +  /Cy  . 


This  approach  requires  no  iteration  and  is  very  efficient;  it  does  involve  several 
additional  matrix  inversions,  however,  since  the  tangent  stiffness  must  be  updated  at 
each  time  step.  For  an  applied  stress  rate  we  may  write  the  Eshelby 
equation  (3-3),  after  substitution,  in  the  form 


(3-9) 


where  AQju  =  Cyu*  -  and  is  the  Kronecker  delta.  Equation  (3-9)  is  of  the 


form 


(3-10) 


Thus, 


®w  ' 


(3-11) 
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With  ki*  in  hand,  one  may  proceed  to  compute  the  rates  of  matrix  stress  and 
composite  strain. 


As  before,  the  matrix  plastic  strain  rate  is  given  by 


ej"  =  F=0  and  ^  0 

=  0  otherwise 


(3-12) 


where  h  is  given  in  Appendix  B  for  two  incremental  plasticity  approaches.  It  is 
important  to  note  that  the  matrix  stress,  plastic  strain,  and  any  associated  state 
variables  (e.g.  backstress)  must  be  updated  to  assign  the  current  hardening  modulus, 
h,  in  the  tangent  stiffness  method. 


3.2  Formulation  for  Comparison  Matrix  Material  in  the  Tangent  Stifbiess  Method 
In  the  general  tangent  stiffness  formulation,  the  comparison  material  may  be 
assumed  to  have  a  compliance  of  the  form 

(CS)-’  ■  -  c^}  *  (C^-'  (^*3) 

or 

(C^)-'  =  (1  »  mc^-'  -  ffCi  (J-M) 

such  that 
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1(0  ’ ojl  i  KO-’oJS  /or  P  ^  0  .  (3-15) 

The  equality  holds  if  /=0  (i.e,,  the  comparison  material  response  is  the  actual  matrix 
material  response),  if  6  =  0,  or  for  purely  elastic  response  of  the  matrix.  Note  that 
(^ijkl  is  the  stiffiiess  of  the  matrix. 

Introduction  of  the  constant  6  is  necessitated  by  the  fact  that  the  response  of 
the  comparison  material  is  not  uniquely  defined  in  the  general  elastic-plastic  loading 
case.  Due  to  the  constraint  associated  with  both  the  fiber  geometry  and  higher  fiber 
stiffness,  the  comparison  matrix  material  generally  has  higher  stress  levels  and, 
consequently,  lower  tangent  stiffness.  It  should  be  emphasized  that  the  concept  of 
comparison  material,  and  of  the  Eshelby  approach  for  that  matter,  is  rooted  in  linear 
elasticity;  hence,  the  incremental  elastic-plastic  case  obviously  must  invoke  the  spirit 
of  the  elastic  approach  but  with  additional  consideration  for  the  actual  degree  of 
constraint  in  the  plastic  regime.  The  parameter  B  affects  the  rate  of  eigenstrain 
accumulation  only  in  the  elastic-plastic  regime.  For  6  =  0,  the  degree  of  constraint 
on  the  matrix  plastic  deformation  is  low  as  is  the  rate  of  eigenstrain  accumulation 
with  composite  strain  accumulation.  As  6  increases,  the  constraint  on  matrix  plastic 
deformation  and,  correspondingly,  the  eigenstrain  accumulation  rate  both  increase. 
This  leads  to  an  increasing  composite  strain  hardening  rate  with  increasing  6. 

The  introduction  of  the  parameter  6  should  therefore  be  viewed  as  a  means 
of  specifying  the  comparison  matrix  material  behavior  in  the  elastic-plastic  regime  in 
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the  incremental  approach.  The  secant  modulus  approach  employed  by  Weng  (1988), 
another  approximate  application  of  the  Eshelby  concept,  does  not  involve  such  an 
assumption  but  is  premised  on  applicability  of  the  deformation  theory  of  plasticity. 
Therefore,  it  is  restricted  to  monotonic,  proportional  loading.  In  general,  it  is 
desirable  to  employ  the  incremental  tangent  stiffness  approach  for  general  loading 
paths.  Moreover,  even  under  proportional  loading  of  the  composite,  the  average 
matrix  plastic  strain  rate  varies  nonproportionally  for  directionally  reinforced 
composites  if  the  applied  loading  is  not  coaxial  with  principal  material  directions. 
This  is  a  first  order  consideration  for  unidirectional  metal-matrix  composites  for  which 
the  plastic  strain  rate  in  virtually  any  off-axis  loading  situation  rotates  away  from  the 
fibers  toward  the  direction  of  least  constraint.  For  such  composite  materials  with 
preferred  directions  of  plastic  flow  which  minimize  constraint,  it  may  be  anticipated 
that  6  will  be  small  since  the  flow  will  seek  these  orientations.  Composites  without 
preferred  reinforcement  directions  (e.g.  particulate  reinforced)  will  require 
significantly  higher  values  of  B  even  though  the  average  level  of  constraint  may  be 
lower  since  there  are  no  preferred  directions  of  plastic  flow  to  minimize  constraint. 
As  a  perhaps  subtle  but  important  consequence  of  the  reinforced  material’s  preferred 
orientation,  one  may  expect  the  possible  existence  of  a  composite  yield  surface  with 
an  associative  flow  rule  to  depend  on  constraint  anisotropy:  the  meaning  of  this  term 
will  become  clear  as  this  discussion  proceeds.  The  term  is  related  to  the  anisotropy 
of  constraint  hardening  of  the  composite,  i.e.  creation  of  internal  stress  associated 


with  inhomogeneity.  At  the  very  least,  such  a  formulation  would  be  highly  complex 
to  model  all  observed  effects  (Wang,  1970).  In  general,  this  is  due  to  the  fact  that 
the  matrix  plastic  strain  rate  rotates  in  a  complicated  manner  relative  to  the 
composite  stress  rate  in  response  to  changing  constraint  during  the  transition  from 
elastic  to  fully  plastic  matrix  response. 

In  this  work,  we  will  apply  the  tangent  stifftiess  approach  to  two  very  different 
classes  of  elastic  reinforcement  in  a  ductile  matrix.  In  the  first  case,  we  will  consider 
a  silica  spherical  particle-reinforced  epoxy  composite  previously  treated  by  Weng 
(1988)  with  the  secant  modulus  approach.  Then,  we  will  consider  unidirectional, 
continuous  fiber-reinforced  boron-aluminum  under  various  loading  conditions. 

In  recognition  of  the  lack  of  constraint  on  matrix  plastic  deformation  for 
fibrous  composites  loaded  off-axis,  we  must  also  introduce  a  provision  for  transition 
to  matrix  plastic  deformation-dominated  behavior.  As  reported  by  Dvorak  and 
Bahei-El-Din  (1987),  the  off-axis  behavior  of  certain  unidirectional  composites 
typified  by  high  fiber  transverse  shear  modulus  is  essentially  dominated  by  matrix 
behavior  in  the  plastic  regime,  nearly  independent  of  fiber  volume  fraction.  For 
composites  with  low  fiber  transverse  shear  modulus,  the  tendency  for  matrix  shear 
localization  and  plastic  strain  rate  rotation  away  from  the  fibers  is  not  as  great. 
Dvorak  and  Bahei-El-Din  have  defined  these  two  types  of  unidirectional  metal-matrix 
composites  as  matrix-dominated  and  fiber-dominated,  respectively. 

While  the  issue  of  defining  the  comparison  material  response  (i.e.  B)  for  the 
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must  be  accomplished  via  a  transition  to  matrix-dominated  deformation  for  off-axis 
loading.  This  transition  behavior  will  depend  both  on  the  plastic  tangent  modulus  of 
the  matrix,  h,  and  the  projection  of  the  direction  of  plastic  straining  in  the  principal 
direction  of  reinforcement  anisotropy  (i.e.  the  fiber  direction),  63.  To  accomplish  a 
decrease  in  the  rate  of  eigenstrain  accumulation  (decreased  misfit  constraint)  as 
plastic  flow  becomes  fully  developed  and  rotates  away  fi’om  the  fiber  direction,  we 
may  introduce  an  effective  fiber  stiffness 


to  be  used  in  the  incremental  BIAS  approach  instead  of  Qj^’  where 


(3-16) 


f  / 

X  =  exp -  -2  1  -  '  . ’  1  -  - 

_[h[  y/2l3  JJ  I  c}^ 


(3-17) 


and  a  and  c  are  the  minimum  and  maximum  diameters,  respectively  of  the  ellipsoidal 
reinforcement.  The  maximum  diameter,  c,  is  assumed  to  be  aligned  with  the  63 
direction.  Constant  h^,  is  a  small  value  of  the  plastic  modulus  h  representative  of  the 
transition  to  asymptotic  (e.g.  past  the  initial  yielding  regime)  matrix  strain  hardening. 
For  the  case  of  long  fibers,  (a/c)-»  0;  (ho/h)'^then  largely  governs  the  transition  to 
matrix-dominated  behavior.  As  (ho/h)’^becomes  much  greater  than  1,  x  0  and  Ciju* 
-»  C,ju^.  This  is  when  matrix-dominated  behavior  is  realized.  Typically,  N  may  be 
expected  to  be  sufficiently  large  (i.e.  N>2)  to  /esult  in  a  continuous  but  rapid 
transition  to  matrix-dominated  behavior.  The  transition  to  matrix-dominated  behavior 
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matrix-dominated  behavior.  As  (ho/h)^  becomes  much  greater  than  1,  x  0  and 
Qjkl^-  This  is  when  matrix-dominated  behavior  is  realized.  Typically,  N  may  be 
expected  to  be  sufficiently  large  (i.e.  N>2)  to  result  in  a  continuous  but  rapid 
transition  to  matrix-dominated  behavior.  The  transition  to  matrix-dominated  behavior 
is  also  affected  by  the  magnitude  of  the  projection  of  the  matrix  plastic  strain  rate  in 
the  fiber  direction,  l^s'nyfijS^j-ejl  =  Injjj.  Clearly,  the  transition  is  delayed  and 
fiber  dominance  persists  until  the  matrix  plastic  strain  rate  rotates  sufficiently  away 
from  the  fiber  direction  (i.e.  |  |  small).  This  condition  is  met  more  rapidly  for 

predominantly  transverse  loadings  than  for  moderate  off-axis  loadings;  hence,  the 
composite  response  remains  fiber-dominated  (Cyy*  ■  Cyy*)  longer  for  loading 
directions  more  nearly  aligned  with  the  fibers. 

^  eigenstrain  rate  diminishes  correspondingly  by  virtue  of 

(3-9),  reflecting  the  loss  of  constraint  on  the  plastic  deformation.  For  fibrous 
composites  (i.e.  low  a/c)  with  a  high  fiber  transverse  shear  modulus,  the  additional 
assumption  Cyu^°  =  Cyy^  (B  =  0)  is  reasonable  since  this  equality  holds  rigorously 
under  both  elastic  and  fully  plastic,  matrix-dominated  conditions. 

For  fibrous  systems  with  low  fiber  transverse  shear  modulus,  it  is  necessary  to 
assign  fi  >  0  since  the  transition  to  matrix-dominated  behavior  is  not  realized  due  to 
the  reduction  of  constraint  anisotropy;  hence,  the  comparison  material  response  is 
slightly  different  from  the  actual  matrix  response.  For  such  cases,  the  rotation  of  the 
matrix  plastic  strain  rate  relative  to  the  fibers  during  loss  of  constraint  anisotropy  and 


interface  shear  localization  is  expected  to  be  less  pronounced  and  associativity  of  a 
flow  rule  based  on  a  composite  yield  surface  is  expected  to  be  more  realistic. 

Note  that  for  spherical  particle  reinforcement,  (a/c)=l  and  %=!  which  leads 
to  =  CjjH*.  Only  B  affects  constraint  hardening  in  the  plastic  regime  in  this  case, 
as  constraint  anisotropy  does  not  enter  as  a  consideration.  The  following  discussion 
will  serve  to  clarify  this  point. 

A  key  idea  in  this  formulation  is  the  concept  of  constraint  hardening.  It  may 
be  obvious  that  for  systems  with  aligned  ellipsoidal  reinforcement,  the  composite  is 
at  most  elastically  transversely  isotropic  if  the  matrix  is  elastically  isotropic.  It  is  well 
known  that  such  fibrous  systems  exhibit  significant  levels  of  "constraint"  on  matrix 
deformation.  A  perhaps  subtle  but  very  important  concept  is  that  of  constraint 
anisotropy  which  has  been  mentioned  previously  in  connection  with  x-  From  a 
formal  standpoint,  it  is  most  instructive  to  view  constraint  in  terms  of  its  evolution 
with  plastic  flow  of  the  matrix.  After  some  level  of  matrix  plastic  straining,  the 
anisotropic  elastic-plastic  state  of  the  composite  may  be  expressed  as  an  isotropic 
tensor  function  of  the  unit  vector  63  in  the  reinforcement  direction  and  backstress  ay® 
associated  with  the  matrix  kinematic  hardening.  The  symmetry  classes  represented 
by  63  and  Oy®  are  transversely  isotropic  and  orthotropic,  respectively.  We  will  refer 
to  the  orientations  of  63  and  ay®  as  structural  and  microstructural,  respectively.  Since 
we  wish  to  examine  the  role  of  plastic  strain  on  constraint  rate  reduction,  we  may 
introduce  the  symmetric  second  rank  tensor  i|ry  defined  in  ey^®  space  by 


(3-18) 
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where  ijTjj  obviously  is  a  vector  in  this  space  normal  to  the  level  surface  cT,jjj“(cyP‘”), 
which  may  be  viewed  as  a  constraint  surface  valid  for  monotonic,  proportional 
loading;  ayj‘“(€,nn'”“)  hydrostatic  constraint  stress  in  the  matrix  which  relaxes 
as  plastic  strains  are  generated.  Clearly,  ijr  jj  is  the  rate  of  increase  of  constraint  stress 
with  plastic  strain.  In  general,  we  may  assign  i|ry  the  form 

This  may  be  expressed  using  the  representation  theorem  in  terms  of  the  irreducible 
functionality  basis  of  the  pair  (^3,  a^®)  (Wang,  1970).  However,  let  us  consider  two 
special  cases  which  are  of  much  relevance  to  both  particle  and  continuous  fiber- 
reinforced  composites. 

First,  for  spherical  particle  reinforcement  and  elastic  isotropy  of  both  phases, 
there  is  no  preferred  structural  orientation  (no  structural  anisotropy)  and  hence 


For  monotonic,  proportional  loading, 


(3-20) 
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provided  the  matrix  has  at  most  only  a  purely  hydrostatic  initial  residual  stress. 
Hence, 
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aearly,  invoking  the  total  derivative  of 


(3-22) 
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The  representation  of  ilfy  is  (Wang,  1970) 


+  ■^3(4»4)^V  ’ 


where 


I,  =  ef  ef  /3  =  ef  e  J-ef 

are  the  two  non-zero  invariants  of  cyP®. 

Hence,  in  general, 

since 

cr  =  A(ii«ri)cr  =  x(4)^r 

for  monotonic,  proportional  loading  and 
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The  constraint  stress  development  is  therefore  independent  of  loading  direction, 
dependent  only  on  the  accumulation  of  plastic  strain  in  the  direction  of  loading.  In 
other  words,  the  constraint  surface  is  isotropic  in  plastic  strain  space. 

For  fibrous  composites,  on  the  other  hand,  structural  anisotropy  dominates; 
and  we  may  consider  the  influence  of  to  be  small  in  comparison  to  that  of  63  with 
respect  to  composite  response.  In  this  case,  we  introduce  the  dependency 


(3-29) 


and  express  the  representation  of  with  respect  to  the  irreducible  functionality  basis 
of  e,,  i.e.. 


'Ifi/  •  <p5i;  + 


(3-30) 


where  <p  >  0,  5  >  0,  and  5  >  <p  such  that  lji|rjjlj  >  0  for  all  vectors,  Ij.  It  is  understood 
that  q>  and  i  may  depend  on  the  invariants  of  0^“  and  the  joint  invariants  of  Oy®  and 
63,  but  the  terms  in  the  representation  based  on  these  generators  have  been  omitted 
as  second  order.  In  this  case,  the  constraint  hardening  rate  for  monotonic  loading 
is  approximated  by 


=  5ieri%  =  5^37. 


(3-31) 


Consequently,  the  constraint  surface  evolves  anisotropically.  Moreover,  the  rate  of 
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constraint  hardening  increases  proportionally  with  nij6j®6j:e3®63=n33,  the  component 
of  the  unit  normal  vector  of  the  plastic  strain  rate  in  the  reinforcement  direction. 

Since 

lifl  ' 

we  may  write  (3-31)  as 

.  m 

—  «  (3J3) 

which  is  the  ratio  of  matrix  constraint  to  projected  deviatoric  stress  rates.  Bear  in 
mind  that  $  depends  on  the  magnitude  of  backstress  as  well.  Note  that  this  analysis, 
based  on  equation  (3-29),  did  not  rely  on  proportionality  of  the  matrix  plastic  strain 
increment  to  achieve  the  result  in  equation  (3-33). 

We  are  now  in  a  position  to  assess  the  form  of  constraint  hardening 
introduced  in  the  incremental  Eshelby  analysis  based  on  these  symmetry 
considerations.  First,  parameter  6  influences  the  comparison  material  response.  For 
the  form  adopted,  for  a  given  value  of  fi,  depends  uniquely  on  and  hence 
h,  the  matrix  plastic  tangent  modulus.  The  constraint  hardening  rate  will  therefore 
depend  on  the  relevant  structural/microstructural  symmetry.  For  spherical  particle 
reinforcement  and  proportionally  applied  tractions. 
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(3-34) 


and  for  fiber  reinforcement  and  proportionally  applied  tractions, 


-1 


(3-35) 


where  the  dependence  of  both  5  and  A’  on  6  is  implicit.  Interestingly,  equation  (3- 
35)  indicates  that  may  become  negative  with  monotonic  plastic  deformation  only 
for  the  fibrous  composite  if  njj  <  0;  physically,  this  condition  is  that  of  a  negative 
matrix  plastic  strain  rate  in  the  fiber  direction.  Clearly,  the  issue  of  distinct  branches 
of  stress-space  composite  yield  criteria  relating  to  matrix-  and  fiber-domination  is 
related  to  the  constraint  anisotropy  inherent  in  the  structural  symmetry  class  of 
transverse  isotropy. 

A  related  issue  is  that  of  %,  which  is  employed  as  a  means  of  transitioning  to 
matrix-dominated  behavior  in  the  incremental  Eshelby  approach,  forcing  the 
eigenstrain  rate  to  vanish  under  fully  plastic  conditions.  This  is  equivalent  in 
principle,  to  the  continuum  slip  model  for  the  matrix-dominated  mode  suggested  by 
Dvorak  and  Bahei-El-Din  (1987)  but  is  also  quite  different  since  a  micromechanical 
framework  is  employed  here  rather  than  construction  of  a  composite  flow  rule.  The 
rate  of  change  of  the  constraint  in  the  Eshelby  approach  is  assumed  to  additionally 
depend  nonlinearly  through  5  on  h  and  njj,  i.e. 
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(3-36) 
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where  again  implicit  dependence  on  6,  h  and  n33  is  expressed.  Note  that  %  employs 
explicit  dependence  on  h  and  njj  to  correlate  off-axis  constraint  weakening  with 
plastic  deformation.  Note  also  that  x  =  0  for  spherical  particle  reinforcement  since 
the  constraint  hardening  is  isotropic  in  this  case.  See  Appendix  F  for  a  more  detailed 
discussion  of  constraint  hardening  in  fiber-reinforced  composites  and  development 
of  the  X  funtion. 


33  Solution  Procedure  for  the  EIAS  Method-Tangent  Stiffness  Formulation 

Per  the  discussion  in  §  3.1  and  3.2,  the  BALNIT  program  was  modified 
accordingly.  The  BALNIT  program  incorporating  the  tangent  stiffiness  method  only 
was  called  TANl.  With  the  exceptions  stated  in  §  3.1,  the  solution  procedure  for 
TANl  remains  the  same  as  BALNIT. 

The  TANl  program  was  modified  to  include  the  6  parameter  and  the  x 
function.  This  program  was  called  CSTAR  since  one  of  its  key  features  is  the 
modification  of  the  fiber  stiffness  (Cy^*)  to  lessen  the  constraining  power.  Again,  the 
solution  procedure  remains  the  same  with  the  exceptions  stated  in  § 
3.2.  Flowcharts  for  these  programs  are  in  Appendix  C. 
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CHAPTER  IV 


DISCUSSION  OF  RESULTS 


4.1  Introduction 

Two  composite  systems  were  investigated  -  boron-aluminum  (B/Al),  a 
unidirectionally  continuous  fiber-reinforced  MMC  and  a  silica  spherical  particle- 
reinforced  epoxy  composite.  The  results  as  applicable  to  B/Al  will  be  discussed  first. 
Discussion  of  the  results  for  the  silica-epoxy  composite  are  in  §  4.7. 

4.2  Application  to  a  Fibrous  Reinforced  Composite 

Due  to  its  ^ood  fiber/matrix  bonding  and  very  periodic  microstructure  with 
virtually  no  presence  of  voids,  B/Al  is  in  many  ways  a  model  material  system  which 
can  be  employed  as  a  standard  to  test  the  accuracy  of  micromechanic  composite 
models  (Pindera  and  Lin,  1989).  The  constituent  material  properties  for  the  B/Al 
system  investigated  are  listed  in  Appendix  E.  These  properties  correspond  to  those 
used  by  Dvorak  and  Bahei-El-Din  in  their  experimental  work  on  B/Al  (1988).  Note 
the  high  transverse  shear  stiffness  of  the  boron  fibers.  The  Eshelby  tensor  for 
circular,  cylindrical  fibers  and  the  nonlinear  kinematic  hardening  rule  utilized  in  this 
investigation  are  shown  in  Appendices  A  and  B,  respectively.  The  actual  uniaxial 


experimental  behavior  of  the  aluminum  matrix  material  as  shown  in  (Pindera  and  Lin, 
1989)  was  used  instead  of  the  "in-situ"  properties  employed  sometimes  in 
development  of  micromechanics  models.  This  in  itself  represents  an  issue  and  a 
possible  advantage  of  the  present  model  over  other  approaches,  as  will  be  discussed 
later.  Uniaxial  experimental  data  on  B/Al,  as  presented  by  Pindera  and  Lin  (1989), 
were  utilized  as  a  basis  for  determining  the  accuracy  of  our  models.  Although  not 
specifically  given,  the  fiber  volume  fraction  of  the  composite  used  in  the  experiments 
appears  to  be  approximately  48%.  We  were  able  to  "back  this  out"  using  the 
constituent  elastic  properties,  the  composite  modulus  for  the  0°  test,  and  a  rule  of 
mixtures  method. 

Uniaxial  histories  were  run  on  the  three  composite  models  (programs 
BALNIT,  TANl,  and  CSTAR)  in  the  0°  and  90°  orientations  as  well  as  three  off-axis 
orientations,  10°,  15°,  and  45°.  These  correspond  to  the  orientations  investigated 
experimentally  by  Pindera  and  Lin  (1989).  The  results  of  the  tests  were  plotted  as 
a^z  versus  where  the  zz-direction  corresponds  to  the  loading  direction.  Note  that 
and  Cj2  represent  the  average  composite  stress  and  strain,  respectively. 

43  Analysis  of  Results  using  the  EIAS  Method  with  Elastic  Constraint 

Figures  4  and  5  show  the  results  using  the  EIAS  Method  with  Elastic 
Constraint.  As  expected,  the  results  are  extraordinarily  stiff.  The  model  does 
reasonably  well  in  the  elastic  regime  since  the  Eshelby  approach  is  based  on  elastic 
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considerations.  However,  once  past  the  point  of  initial  yielding,  continued 
accumulation  of  matrix  plastic  strain,  which  is  treated  as  a  negative  eigenstrain,  causes 
increased  constraint  and  an  overly  stiff  response.  These  results  confirm  what  has 
been  stated  consistently  in  the  literature,  that  unmodified  self-consistent  approaches 
are  inadequate  for  modeling  elastic-plastic  behavior  of  composites,  particularly  for  the 
case  of  the  ductile  matrix. 

4.4  Analysis  of  EIAS  Method-Tangent  Stiffiiess  Formulation 

Figures  6  and  7  show  the  results  of  the  tests  using  the  EIAS  Method-Tangent 
Stiffness  Formulation  with  B=0  and  without  the  %  function,  i.e.  always. 

With  6=0,  the  comparison  material  is  the  same  as  the  matrix  material.  The  model 
does  reasonably  well  in  approximating  initial  yielding  behavior  but  does  not  capture 
the  transition  to  asymptotic  plastic  flow.  Without  any  correction  for  unwarranted 
eigenstrain  accumulation,  we  see  the  composite  stiffen  in  the  plastic  regime  following 
initial  yielding.  These  results  demonstrate,  however,  that  the  tangent  stiffiiess  method 
has  improved  the  model’s  ability  to  approximate  initial  yielding  (i.e.  the  position  and 
shape  of  the  "knee"  in  the  curve).  However,  eigenstrain  accumulation  or,  conversely, 
excessive  constraint  hardening,  still  presents  a  problem.  Extensive  analysis  was 
conducted  to  determine  the  effect  of  various  values  of  6  on  the  response  and  to 
determine  if  6  alone  would  reduce  the  effects  of  excessive  eigenstrain  accumulation 
and  bring  the  response  more  in  line  with  the  experimental  data.  These  efforts  proved 
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futile,  as  would  be  expected  based  on  the  discussion  of  constraint  anisotropy  in  the 
previous  section.  A  particular  value  of  B  would  suffice  only  for  a  particular  loading 
direction.  It  was  evident  that  the  %  function  is  necessary  to  model  the  change  of 
constraint  hardening  for  all  orientations.  The  6  factor  is  a  comparison  material 
parameter  which  should  be  unique  to  a  particular  composite  system  and  should  not 
change  with  each  loading  orientation  for  even  continuous  fiber-reinforced  composites. 

4.5  Analysis  of  Rotation  of  the  Plastic  Strain  Rate  Rector 

Based  on  the  literature,  the  inelastic  response  of  the  composite  should 
asymptotically  approach  (or  closely  approximate)  that  of  the  pure  matrix  material  for 
certain  ranges  of  off-axis  uniaxial  loading  directions  relative  to  the  fiber  direction. 
This  motivated  an  investigation  of  the  rotation  of  the  plastic  strain  rate  vector,  CjjP"*, 
as  a  means  of  studying  the  transition  from  fiber-dominated  to  this  so-called  matrix- 
dominated  behavior. 

To  help  illustrate  this  point,  we  employed  the  tangent  stiffness  method  with 
6=0  (no  X  function)  and  plotted  the  matrix  plastic  strain  in  the  fiber  direction,  CajP™, 
versus  the  matrix  plastic  strain  in  the  transverse  direction,  622*’“  for  several  different 
monotonic,  uniaxial  loading  orientations  of  interest.  The  directions  33  and  22  refer 
to  the  fixed  material  coordinate  system,  33  being  the  direction  of  the  aligned  fibers. 
These  results  are  shown  in  Figure  8.  It  is  important  to  note  that  the  tangent  stiffness 
method  accurately  modeled  the  initial  yielding  regime.  Therefore,  it  is  accurate  up 
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to  approximately  0.2%  strain.  Therefore,  the  results  shown  in  Figure  8  are  also 
accurate  for  this  same  strain  range.  For  the  0°  orientation,  relatively  significant  levels 
of  positive  633*”“  are  generated  while  only  limited  levels  of  negative  622'”“  are 
observed.  This  is  in  accordance  with  plastic  incompressibility  of  the  matrix.  No  CyP*" 
rotation  takes  place  in  the  0®  direction.  In  the  10°  orientation,  the  response  is  initially 
very  similar  to  that  of  the  0°  orientation;  however,  rotation  of  CyP®  causes  generation 
of  positive  €22*’“,  and  we  see  the  curve  rotate  towards  positive  622*’“-  In  the  15° 
orientation,  we  observe  behavior  similar  to  the  10°  orientation.  Due  to  faster  rotation 
of  CyP“,  the  curve  turns  more  rapidly  to  the  right,  again  demonstrating  the  increasing 
dominance  (or  positive  magnitude)  of  C22''"”.  The  behavior  of  the  45°  and  90°  loading 
angles  are  essentially  matrix-dominated  from  the  outset.  The  negative  values  of  C33P® 
observed  in  the  90°  direction  are  realized  due  Poisson’s  effect.  Based  on  these 
results,  we  recognize  the  tendency  of  the  matrix  plastic  strain  rate  vector  to  rotate 
and  align  itself  with  the  22-direction.  Since  predominance  of  822**“  in  the  response 
is  indicative  of  matrix-dominated  behavior,  we  conclude  that  the  transition  to  matrix- 
dominated  behavior  is  related  to  the  rotation  of  CyP“.  Only  a  micromechanical  model 
can  accurately  predict  this  effect.  It  is  very  important  to  note  that  this  effect  is 
predicted  without  introduction  of  the  %  function;  x  merely  serves  to  affect  the 
eigenstrain  rate  past  the  initial  yielding  regime.  Also,  since  the  matrix  plastic  strain 
rate  is  nonproportional  even  though  the  tractions  are  applied  proportionally,  this 
casts  grave  doubts  on  the  possibility  of  constructing  a  simple,  associative  flow  rule 
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plasticity  framework  based  on  a  composite  yield  surface,  at  least  for  composites  which 
exhibit  this  transition  to  matrix-dominated  behavior, 

4.6  Analysis  of  the  EIAS  Method-Tangent  Stiffness  Formulation  with  0  Parameter 
and  X  Function 

Figures  9  and  10  show  the  results  of  the  tangent  stifrness  method  including  the 
X  function.  6=0  was  selected  on  the  basis  that  the  B/Al  composite  is  expected  to 
exhibit  a  matrix-dominated  regime  where  the  incremental  composite  response  is 

1 

essentially  governed  by  that  of  the  matrix.  In  the  %  function,  the  parameters  ho=  180 
ksi  and  N=4  were  selected  to  best  fit  the  uniaxial  test  in  the  45°  orientation.  Good 
*  results  are  realized  in  all  orientations  as  the  model  now  approximates  a  continuous 

transition  to  matrix-dominated,  plastic  deformation  for  off-axis  loading.  The  x 
function  controls  the  rate  of  constraint  hardening;  therefore,  we  do  not  see  increased 

I 

constraint  and  the  associated  stiffening  caused  by  eigenstrain  accumulation  as  was  the 
case  with  the  model  without  the  x  function.  Comparing  the  10°  response  to  the  45° 
I  response,  the  reduction  of  constraint  is  more  gradual  in  the  10°  orientation;  the 

model  "recognizes"  that  the  loading  is  more  nearly  aligned  with  the  fiber  direction 
and  hence  the  response  will  remain  fiber-dominated  longer  (i.e.  it  takes  longer  to 

I 

realize  matrix-dominated  behavior). 

It  is  important  to  point  out  that  use  of  an  incremental  plasticity  theory  enabled 
^  us  to  conduct  the  rotation  analysis  of  6^*””.  The  secant  modulus  method  (Weng, 
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1988),  in  contrast,  cannot  support  such  studies  because  a  backstress  computation, 
inherent  to  an  accurate  incremental  theory,  is  required  to  analyze  ejjP'"  rotation  with 
some  degree  of  accuracy.  Also,  it  is  anticipated  that  the  constraint  anisotropy  issue 
will  not  be  as  pronounced  for  fibrous  composites  with  low  transverse  elastic  shear 
modulus;  in  this  case,  a  combination  of  6>0  and  the  %  function  is  most  likely 
necessary. 

The  goal  of  concentric  cylinder  models  which  incorporate  interface  phases 
(Mikata  and  Taya,  1985)  as  well  as  the  recently  proposed  "smear"  model  (Luo  and 
Weng,  1989)  both  seek  to  describe  the  interface  shear  localization  and  loss  of 
constraint  by  modifying  the  transformation  strains.  In  contrast,  our  approach  is  to 
introduce  an  effective  fiber  stiflbiess;  nonetheless,  the  goal  is  similar. 

4.7  Application  to  Particle-Reinforced  Composites 

Next,  we  apply  the  tangent  stiffness  model  to  the  elastic-plastic  behavior  of  a 
silica  particle-reinforced  epoxy  composite.  The  silica  is  elastic,  and  the  epoxy  is 
elastic-plastic.  Both  phases  are  assumed  to  be  isotropic.  The  material  properties  are 
from  Weng  (1988)  and  are  listed  in  Appendix  E.  The  Eshelby  tensor  for  spherical 
particle-reinforcement  and  the  isotropic  hardening  plasticity  theory  employed  to 
model  the  epoxy  matrix  behavior  are  given  in  Appendices  A  and  B,  respectively. 

As  stated  in  section  32,  there  is  no  structural  anisotropy  with  spherical 
particle-reinforced  composites.  The  fiber  diameter  ratio  (a/c)  =  1;  therefore,  %  = 
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1,  and  Cjjy*  =  Cjjy*.  Hence,  only  6  affects  constraint  hardening.  The  BIAS  Method- 
Tangent  Stiffness  Formulation  with  the  6  parameter  for  calculation  of  comparison 
matrix  material  response  was  used  to  calculate  the  uniaxial  response  of  the  composite 
with  the  same  volume  fractions  as  in  Weng’s  work  (1988). 

Experimental  data  (Weng,  1988)  were  digitized  and  plotted  with  the  results 
of  the  tangent  stiffness  model  in  Figures  11-13.  The  composite  stress  versus  plastic 
strain  is  considered  in  these  plots  and  is  of  course,  orientation  independent  for  each 
volume  fraction  of  reinforcement. 

The  composite  response  was  calculated  with  6=0  and  6=8.  The  plot  for  6=0, 
Figure  11,  represents  behavior  of  the  composite  when  the  comparison  matrix  material 
response  is  the  same  as  that  of  the  actual  matrix  material.  The  results  of  the  tangent 
stiffness  model  best  fit  Weng’s  secant  modulus  approach  when  6=8.  These  results 
are  plotted  in  Figure  12,  while  the  experimental  results  are  shown  in  Figure  13. 

Several  points  are  noteworthy.  First,  the  results  agree  relatively  well  with  the 
experimental  data.  Second,  the  results  are  very  similar  to  those  of  Weng  who 
employed  the  EIAS  method  with  a  secant  modulus  based  on  the  deformation  theory 
of  plasticity.  It  therefore  appears  that  the  tangent  stiffness  method  is  a  reasonable 
extension  of  the  secant  modulus  concept  for  the  solution  of  incremental  elastic-plastic 
problems  involving  complex  cyclic  loading  paths.  Refer  to  Appendix  D  for  a  more 
complete  comparison  of  the  two  methods. 

Due  to  the  linear  combination  of  6  and  f  in  (3-14),  the  response  stiffens  for 
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a  particular  B  as  the  volume  fraction  /  increases.  An  increase  in  /  leads  to  increased 
constraint  and,  consequently,  a  much  stiffer  response  in  the  yielding  regime.  This  fact 
^  is  borne  out  by  Figures  11  and  12.  The  stress  versus  plastic  strain  curves  with  6=0 

are  much  flatter  than  those  with  6=8  for  the  same  volume  fraction  of  reinforcement. 
Finally,  our  results  with  the  tangent  stiffness  method  with  6=8  are  stiffer  than 
^  Weng’s  results,  especially  at  lower  volume  fractions.  One  notes,  however,  that  the 

actual  experimental  results  are  also  stiffer  than  Weng’s  theoretical  results  (Weng, 
^  1988). 

The  material  response  calculated  using  the  BIAS  model  with  elastic  constraint, 
as  shown  by  Weng  (1988),  is  significantly  stiffer  than  experimentally  observed.  This 
^  concurs  with  the  results  of  our  earlier  calculations  on  B-AJ  fibrous  composites. 

I 

I 

I 
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CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

5.1  Summary  of  Composite  Material  Model 

As  a  result  of  our  theoretical  study  of  the  inelastic  behavior  of  unidirectional, 
continuous  fiber-reinforced  metal-matrix  composites,  we  developed  a  practical, 
computationally  efficient  micromechanical  model  which  satisfactorily  approximates 
experimental  inelastic  behavior  of  certain  composites.  This  model  is  limited  to  two- 
phase  composite  systems  where  the  matrix  is  the  ductile  phase  but  can  be  modified 
for  any  aspect  ratio  of  aligned  reinforcement  by  employing  the  appropriate  Eshelby 
tensor.  Material  properties  of  the  constituent  phases  must  be  input  to  the  program. 
Detailed  information  regarding  the  stress  and  strain  states  of  the  composite  as  well 
as  that  of  the  individual  phases  can  be  obtained  from  the  model. 

The  model  employs  a  sophisticated  nonlinear  kinematic  hardening  rule  to 
model  the  plastic  flow  of  the  matrix.  Use  of  this  incremental  plasticity  theory 
necessitates  introduction  of  a  tangent  stifiness  method  instead  of  a  secant  stiffness 
method,  which  is  restricted  to  deformation  plasticity  theory.  Furthermore,  use  of  an 
incremental  theory  facilitates  consideration  of  nonproportional  and/or  cyclic  loading 


histories. 


I 


The  model  employs  two  additional  parameters,  6  and  %,  to  account  for  the 
constraint  effects  associated  with  eigenstrain  accumulation  and  anisotropy  due  to  fiber 
reinforcement.  Although  no  exact  function  for  the  B  parameter  has  been  determined, 
intuitively  it  is  u  function  of  reinforcement  volume  fraction,  fiber  geometry,  and 
elastic  properties  of  the  reinforcement.  The  B  parameter  is  most  singularly 
applicable  to  composite  systems  with  little  or  no  reinforcement  or  structural 
anisotropy  while  the  %  function  addresses  cases  of  structural  anisotropy,  e.g. 
unidirectional,  continuous  fiber-reinforced  systems. 

5.2  Significant  Findings 

Due  to  the  heterogeneous  character  of  composites,  the  behavior  of  one 
constituent  is  not  mutually  exclusive  of  that  of  the  others.  This  interaction  of  the 
phases  renders  analysis  of  the  homogeneous  aggregate  very  difficult,  particularly  in 
the  plastic  regime.  In  this  thesis,  we  recognize  the  deficiencies  of  the  EIAS  method 
with  elastic  constraint  and  pursue  sound  methods  of  modifying  the  comparison 
material  to  collectively,  and  perhaps  indirectly,  take  into  account  the  first  order 
effects.  Our  ultimate  goal  in  this  process  was  to  effectively  model  the  inelastic 
behavior  of  unidirectional,  continuous  fiber-reinforced  composites. 

In  parallel  to  Weng’s  secant  modulus  method  (1988),  we  introduced  the 
tangent  stiffness  method.  Our  use  of  an  incremental  plasticity  theory  enabled  us  to 
instantaneously  compute  the  matrix  plastic  hardening  modulus  which  was  used  to 
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compute  the  "tangent  stiffness  modulus."  This  method  produced  equal,  if  not 
better,  results  than  Weng’s  secant  modulus  method  based  on  a  comparison  for  a 
spherical  particle-reinforced  composite  system,  but  failed  to  effectively  approximate 
the  experimental  behavior  for  fibrous  systems. 

Recognizing  that  the  overly  stiff  response  was  due  to  the  overconstraining 
effects  unique  to  fiber  reinforcement,  we  focused  on  the  realization  that  experimental 
results  evidenced  an  eventual  but  complete  transition  to  matrix-dominated  behavior 
for  all  off-axis  loadings,  at  least  for  the  B/Al  composite  considered.  From  the  analysis 
in  §  4.5,  we  identified  the  rotation  of  the  matrix  plastic  strain  rate  vector  as  the 
primary  cause  for  the  transition  to  matrix-dominated  behavior.  Our  %  function 
employs  this  rotation  concept  as  a  means  to  reduce  the  effect  of  fiber  constraint 
anisotropy  and  represents  a  significant  development  above  and  beyond  any  methods 
currently  in  existence  for  modeling  inelastic  behavior  of  fibrous  composites.  By 
incorporating  the  x  function  into  the  EIAS  method-tangent  stiffness  formulation,  we 
successfully  approximated  the  behavior  of  a  fibrous  system. 

We  also  introduced  the  6  parameter  as  a  means  of  accounting  for  the  isotropic 
constraint  effects  in  composites.  This  parameter  is  more  relevant  to  composite 
systems  with  little  or  no  reinforcement  anisotropy  (e.g.  spherical  particle-reinforced 
composites)  instead  of  systems  reinforced  by  unidirectional,  continuous  fibers.  As 
stated  in  §  3.2,  fibrous  composites  may  have  distinct  preferred  directions  of  matrix 
plastic  flow  due  to  their  aligned  reinforcement  whereas  there  are  no  such  clearly 
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defined  directions  of  flow  in  the  spherical  particle-reinforced  composite.  Additionally, 
the  fibrous  system  with  high  fiber  shear  stiffness  will  tend  to  reduce  constraint 
hardening  rate  through  the  phenomenon  mentioned  above,  matrix  plastic  strain  rate 
rotation,  as  well  as  interface  shear  localization.  The  particulate  system  cannot  relieve 
constraint  preferentially;  therefore,  it  experiences  additional  constraint  in  the 
asymptotic  plastic  flow  regime.  Hence,  the  6  parameter  is  required  for  the  spherical 
particle-reinforced  system,  but  it  is  not  needed  with  the  fibrous  system  considered. 

Perhaps  one  of  our  most  significant  findings  is  that  we  have  substantiated  the 
experimental  observations  of  B/Al  behavior  made  by  Dvorak  (1987)  which  formed 
the  basis  of  their  bimodal  plasticity  theory  of  fibrous  composites  (Dvorak  and  Bahei- 
El-Din,  1987).  Due  to  the  plastic  strain  rate  rotation  exacerbated  by  the  high  boron 
fiber  transverse  shear  stiffness,  we  support  Dvorak  and  Bahei-El-Din’s  findings  of 
non-normality  of  the  plastic  strain  rate  vector  and  "flats"  in  the  composite  yield 
surface.  However,  we  feel  that  it  is  inadvisable  to  attempt  to  model  a  composite  such 
as  B/Al,  which  so  clearly  demonstrates  both  modes  of  deformation,  with  a  normality 
flow  rule  based  on  a  composite  aggregate  yield  surface.  This  type  of  approach  might 
be  more  successful  when  applied  to  a  system  with  more  compliant  fibers  in  transverse 
shear.  We  find  it  particularly  interesting  that  we  have  recovered  Dvorak  and  Bahei- 
El-Din’s  anomalous  results  through  an  entirely  different  approach  based  on 
micromechanics. 

Finally,  we  find  it  quite  interesting  that  we  were  able  to  achieve  good  results 
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without  consideration  of  the  so-called  "in-situ"  behavior  of  the  matrix.  We  utilized 
the  actual  uniaxial  experimental  behavior  of  the  aluminum  matrix  material  as  shown 
in  Pindera  and  Lin  (1989)  to  determine  the  material  parameters  required  to  model 
matrix  flow  instead  of  the  "in-situ"  properties  used  by  so  many  in  the  development 
of  their  micromechanics  models  (Pindera  and  Lin,  1989).  This  in  itself  represents  a 
significant  departure  from  contemporary  approaches  and  may  be  viewed  as  an 
advantage  of  our  model  over  these  approaches,  unless  the  "in-situ"  behavior  is  (i) 
warranted  based  on  different  dislocation  structure  and/or  (ii)  determined  to  coincide 
with  that  actually  employed  in  other  models.  Otherwise,  the  assumed  matrix  response 
in  other  models  may  be  viewed  as  a  kind  of  "first  order  fitting  factor."  As  an 
example,  Pindera  and  Lin  (1989)  adjusted  the  actual  material  properties  to  the  so- 
called  "in-situ"  properties  to  allow  for  the  effects  of  fabrication  on  the  matrix  material 
and  so  that  their  analytical  results  correlated  with  the  experimental  data.  The  matrix 
response  is  'Tracked  out"  by  fitting  the  composite  response.  This  'Tsacking  out" 
process  is  not  required  with  our  model. 

5J  Recommendations 

As  stated  before,  modeling  of  composite  behavior  is  very  complex  because 
there  are  so  many  variables  to  take  into  consideration.  For  instance,  shear- 
localization  at  the  fiber-matrix  interface  is  universally  recognized  as  having  a 
significant  effect  on  the  strength  of  fibrous  composites,  especially  systems  such  as 
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B/AI  which  has  a  high  fiber  transverse  shear  modulus.  However,  we  achieved 
acceptable  results  without  specifically  addressing  this  in  our  model.  Perhaps  we 
indirectly  accounted  for  it  in  our  introduction  of  %  for  structural  constraint  anisotropy. 
Application  of  our  model  to  systems  reinforced  by  compliant  fibers  is  outside  the 
scope  of  this  thesis.  Therefore,  we  have  several  recommendations  for  future  research 
to  further  test  the  validity  of  our  model. 

First,  the  EIAS  model-tangent  stiffiiess  method  with  the  %  function  should  be 
applied  to  a  system  reinforced  by  compliant  fibers  (eg.  graphite)  to  determine  how 
well  it  replicates  experimental  data.  Certainly,  we  expect  the  %  function  to  play  a 
different  role  for  the  compliant  system. 

Secondly,  it  would  be  interesting  to  conduct  a  similar  plastic  strain  rate 
rotation  analysis  on  this  compliant  system  to  see  if  it  demonstrates  both  matrix-  and 
fiber-dominated  deformation  modes.  According  to  Dvorak  and  Bahei-El-Din,  a 
system  such  as  graphite-aluminum  should  deform  only  in  the  fiber-dominated  mode 
(1987).  If  this  is  the  case,  plots  of  €33*””  versus  622’’“  for  uniaxial,  monotonic  loading 
in  off-axis  orientations  should  retain  a  higher  degree  of  proportionality  than  that  in 
Figure  8. 

Finally,  we  stated  that  we  utilized  an  incremental  theory  so  that 
nonproportional  and/or  cyclic  loading  histories  could  be  applied  to  our  model.  We 
did  not  report  any  such  histories,  but  it  would  be  interesting  to  see  how  our  model 
correlates  with  the  data  recorded  by  Dvorak  and  Bahei-El-Din  for  nonproportional 
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loading  experiments  conducted  on  B/Al  (1988). 
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Infinite  material  domain  D  containing  an  ellipsoidal  inhomogeneity,  n, 
lied  stress  0^°  (a)  actual  problem,  (b)  equivalent  inclusion  problem 
Arsenault,  1989). 
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Figure  4.  Results  of  the  EIAS  method  with  elastic  constraint  for  B/Al  -  0°  and  90° 
orientations. 


Figure  5.  Results  ot  the  EIAS  method  with  elastic  constraint  for  B/Al  -  10°,  15°,  and 
45°  orientations. 
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Figure  6.  Results  of  the  BIAS  method-tangent  stiffness  formlu'>ition  without 
parameter  or  x  function  for  B/AJ  -  0“  and  90®  orientations. 
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Figure  7.  Results  of  the  EIAS  method-tangent  stiffness  formulation  without  fi 
parameter  or  x  function  for  B/Al  -  10®,  15®,  and  45®  orientations. 
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Figure  8.  Analysis  of  matrix  plastic  strain  rate  rotation  for  0®,  10®,  15®,  45®,  and 
orientations. 
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Figure  9.  Results  of  the  BIAS  method-tangent  stiffness  formulation  with  6=0  and 
the  X  function  for  B/AI  -  0®  and  90®  orientations. 
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Figure  10.  Results  of  the  BIAS  method-tangent  stiffness  formualtion  with  fl=0  and 
the  X  function  for  B/Al  -  10®,  15®,  and  45®  orientations. 


Results  of  the  EIAS  method-tangent  stiffness  formulation 
spherical  particle-reinforced  composite. 


APPENDIX  A 

ESHELBY’S  TENSOR  OF  ELASTICITY,  Syy 

A.1  Introduction 

In  the  context  of  the  equivalent  inclusion  problem,  the  Eshelby  tensor  of 
elasticity,  8^^,  relates  the  eigenstrain,  Cy*,  to  the  perturbance  strain,  In  general, 
Sjjy  is  a  nonsymmetric  tensor  calculated  based  on  the  geometry  of  the  composite 
reinforcement.  Mura  (1987)  derives  the  Eshelby  tensor  in  great  detail  for  both 
isotropic  and  anisotropic  inclusions.  Since  the  composite  systems  analyzed  in  this 
thesis  had  isotropic  inclusions,  we  show  here  the  essentials  to  understand  Mura’s 
derivation  of  the  Eshelby  tensor  for  isotropic  inclusions  and  give  the  equations  of  the 
Eshelby  tensor  for  spherical  and  continuous  fiber-reinforced  composites.  The 
problem  and  associated  field  equations  are  as  set  forth  in  §  2.1  and  §  2.2  of  the  main 
thesis. 

A.2  General  EIxpressions  of  Elastic  Fields  for  Given  Eigenstrain  Distributions 

Eigenstrain  is  introduced  as  a  description  of  misfit  due  to  periodic 
inhomogeneity  typical  of  reinforced  composites.  The  corresponding  eigenstresses  are 
self-equilibrated. 


Following  Mura  (1987),  the  fundamental  linear  elastic  equations  to  be  solved 


for  given  eigenstrain  Cy*  are 

where  Uj  is  the  displacement  at  an  arbitrary  point  x(xiPC2,X3).  The  body  is  assumed 
to  be  free  of  any  external  surface  traction  or  body  force  such  that  ayj=0  everywhere 
within  the  body  and  aynj=0  on  the  boundary.  We  will  consider  the  body  to  be 
infinitely  extended  such  that  Uy  -►  0  as  Xj  «  reproduces  the  traction  free  condition. 
In  the  case  of  periodic  solutions,  suppose  ey*(x)  is  given  in  the  form  of  a  single  wave 
of  amplitude  Cy*(5),  such  as, 

cj(*)  *  (^-2) 

where  5  is  the  wave  vector  corresponding  to  the  given  period  c '  the  distribution,  and 


i  =  and 

since  5  =  5iSi-  The  solution  of  (.4-1)  corresponding  to  this  distribution  is  obviously  in 
the  form  of  a  single  wave  of  the  same  period,  i.e. 

=  5X5)exp(i5*Jf).  (A-4) 

Substituting  (A-2)  and  (A-4)  into  (A-1),  we  have 


~  (A  5) 

which  represents  three  equations  for  determining  the  three  unknowns  Uj  for  a  given 
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Using  the  notation 


«»(0  =  V/. 

the  displacement  amplitude  vector  Qj  is  simply  obtained  as 

where  Njj  are  the  cofactors  of  the  matrix  Kij(0  and  0(5)  is  the  determinant  of  KyCO- 
Substituting  (A-6)  into  (A-4),  we  have 

The  corresponding  equations  for  strain  and  stress  are 


-  C,*^C^r^(5)5,5A^(Oexp(«5.*) 
-  ««(*)}• 
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A3  Method  of  Fourier  Integrals 

A  more  sophisticated  distribution  of  eigenstrain  than  the  single  wave  in 
equation  (A-2)  may  be  introduced.  Of  course,  any  distribution  may  be  represented 
by  the  single  wave  solution  in  the  Fourier  series  sense.  If  is  given  in  the  Fourier 
integral  form  for  an  infinite  domain,  namely, 

m 

—m 

where 


J'/O  =  (2nr’  /  e;<i)(a[p(-i5.x)rfr,  (A-10) 

then  the  displacement,  stress,  and  strain  can  be  expressed  as 


and 


«Xx)  =  -C")-’;!-  /  /  AjCSW’-'CO 

-m  -«• 

xexp{i^»(jr  -  x')}d^dx', 

=  (2n)-’/  /  •„(*') 

X  exp{i5*(x  -  x')]d^dx'. 


(A-lla) 
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-m 

xexp{i5*(x  -  x')}d^dx'  -  c^x)| 


for  general  eigenstrain  distribution  Cij*(x). 


A.4  The  Method  of  Green’s  Functions 


When  Green’s  functions  Gy(x-x’)  are  defined  as 


G^<x-x')  =  (2izr^  f  N^a)D-H^)ejpm»(x-x')}d^, 


the  displacement  Uj(x)  in  (A-11)  can  be  rewritten  as 


«.(x)  =  -  f  Cj^tl,(x')GJx  -  x')dx', 


where  Gy  ,(x-x’)  =  3/6!x,{Gy(x-x’)  =  -5/5x,’{Gjj(x-x’). 

The  corresponding  expressions  for  the  strain  and  stress  become 


•• 
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iiu 


(A.15) 


In  1963,  Mura  (1987)  rewrote  aj:(x)  in  (A-15)  in  the  form 


«(,<*)  =  Ctftf  /  (A-16) 


where  e^^  are  permutation  symbols.  Since  esihe,„,,  =  5si<^tn‘'^sn^ti»  (A-16)  becomes 


(A-17) 


Following  Mura  (1987),  it  is  shown  that 

(A-IS) 

where  5(x-x’)  is  Dirac’s  delta  function  having  the  property 


/  tj,x')b{x-x')dx'  =  c;;^x). 

■^m 

It  is  seen  from  (A-18)  that  the  Green’s  function  Gpfc(x-x’)  is  the  displacement 
component  in  the  Xp-direction  at  point  x  when  a  unit  body  force  in  the  x^-direction 
is  applied  at  point  x’  in  the  infinitely  extended  material.  By  this  definition  of  the 
Green’s  function,  we  can  derive  (A-13)  directly  from  (A-1).  As  mentioned  previously, 
the  displacement  Uj  in  (A-1)  can  be  considered  as  a  displacement  caused  by  the  body 
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force  ,*  applied  in  the  Xj-direction.  Since  Gjj(x-x’)  is  the  solution  for  a  unit 

body  force  applied  in  the  Xj-direction,  the  solution  for  the  present  problem  is  the 
product  of  Gjj  and  the  body  force  namely, 


u.(x)  =  -  f  G^ix-x')Cj^tU;c')dx'. 


(A-19) 


Integrating  by  parts  and  assuming  that  the  boundary  terms  vanish,  we  have. 


“iW'  /  (A-2«) 


Expression  (A-19)  is  preferable  to  (A-20).  When  is  constant  in  n  and 
zero  in  D-n  (as  per  the  Eshelby  solution),  it  can  be  seen  that  the  integrand  in  (A-20) 
vanishes  except  on  the  boundary  of  n  since  Cnw**® 


A.5  Explicit  Forms  of  the  Green’s  Function  for  Isotropic  Inclusions 

An  inclusion,  containing  eigenstrain  is  isotropic  if  it  has  the  same  elastic 
moduli  as  its  surrounding  domain.  For  isotropic  materials,  the  integrands  in 
equations  (A-11)  are 


(A-21) 


D(0  =  iiHX  ^  2p)5«, 

^  2p)6^,.52  -  ik  + 

where  and  A,  and  n  are  Lame’s  elastic  constants  (n  is  the  shear  modulus). 
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If  we  substitute  the  expressions  in  (A-21)  into  the  integral  expression  for  the 
Green’s  Function  (A-12)  and  integrate,  we  obtain 


°76n^a^v)|x|t<^  ■  • 


(A.22) 


where  |x|  or 


G^x-x^  = 


1  \ 


^\x-x' 


(A-23) 


4Ttp  (jc-x^l  1671  p(l  -  v)  ^x^^XJ 
where  v  is  Poisson’s  ratio  and  |  x-x’  j  =(Xi-Xi’)(Xj-Xj’).  Equation  (A-22)  or  (A>23)  is  the 
form  of  the  Green’s  function  we  need  to  obtain  explicit  expressions  for  the  Eshelby 
tensor  based  on  the  geometry  of  the  reinforcement. 

From  (A-20),  we  have  the  following  expression  for  the  displacement  in  terms 
of  the  Green’s  function: 


(A-24) 


where  n  is  an  ellipsoidal  inclusion  given  by  the  equation 


2  2  3 

X,  X2  Xj 

—  +  —  +  —  si, 
2  2  2 

fll  <22  «3 


where  a^,  a2>  ^3  ellipsoid. 

Gij(x-x’)  is  from  (A-22). 
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After  some  manipulation,  equation  (A-24)  can  be  rewritten  as 


u,(x)  =  - — — f  gM(n) — - 

‘  871(1  -  v)JQ  I-/ 


dx' 


x'-x\ 


(A-25) 


where 


%(»)  =  (1  -  2v)(6^4  * 

The  vector  n  is  the  unit  vector  (x’-x)/|  x*-x  j . 

When  point  x  is  located  inside  the  inclusion,  the  integral  in  (A-25)  can  be 
explicitly  performed.  The  volume  element  dx’  in  (A-25)  can  be  written  as 
dx’=drdS=drr^do),  where  r  =  |  x’-x  |  and  dw  is  a  surface  element  of  a  unit  sphere  E 
centered  at  point  x.  Upon  integration  with  respect  to  r,  we  have 


(A-27) 


87c(1  - 

Following  Mura  (1987),  the  integral  (A-27)  may  be  transformed  to 


u^ix)  = 


x^e 


87t(l  -  g 


Jt 


^da>. 


(A-28) 


where 


and 
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3  2 
<-i  af 


From  Cij(x)=y2(Uij(x)+Uj  j(x)),  the  strain  components  become 


e  <x) - ^ - f 

^  16x(l  -  v)Ji:  g 

The  integral  in  (A-29)  is  independent  of  x.  Therefore,  we  have  verified  Eshelby’s  key 
result  in  that  the  strain  (and  therefore  the  stress)  is  uniform  inside  the  inclusion. 

It  is  convenient  to  write  (A-29)  as 


(A-30) 


where  is  called  the  Eshelby  tensor;  Sjjy  takes  on  the  general  form 


~  16x(l  - 


1 _ 

[  _  v)Jj:  g 


(A-31) 


for  an  isotropic  material  with  an  inclusion.  Note  that  the  Lame’s  constants 
correspond  to  the  matrix  since  the  isotropic  inclusion  is  assumed  to  have  the  same 
elastic  properties  as  the  matrix  but  with  an  eigenstrain  introduced  in  n.  The  Eshelby 
tensor  is  not  symmetric  in  general. 


A.6  Eshelby’s  Tensor  for  Spherical  Inclusions 

In  our  analysis  of  spherical  particle-reinforced  composites,  we  utilize  the  form 
of  the  Eshelby  tenser  for  spherical  inclusions.  When  the  inclusions  are  spherical  in 
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shape,  ai=a2=a3,  and  (A-31)  simplifies  dramatically.  The  equations  for  the  non-zero 
components  of  the  Eshelby  tensor  are: 


“  *^2222  “  ^3333 


7  -  5v 
15(1  -  V) 


^1122  “  ^2233  ~  ^3311  ~  ^1133  "  ^2211  “  ^3322 


5v  -  1 
15(1  -  V)  ’ 


(A.32) 


^1212  “  ^2323  “  ^3131 


4  -5v 
15(1  -  V) 


A.7  Eshelby’s  Tensor  for  Continuous  Fibers 

If  the  inclusions  are  elliptical  cylinders  (i.e.,  -»  «)  as  is  the  case  with 

continuous  fibers,  then  the  non-zero  components  of  the  Eshelby  tensor  are: 


’nil 


'2222 


’1122 


*2233 


'1133 


'2211 


1 

o^  +  2a^a^ 

2(1  -  V) 

(Oj  +  Oj)^ 

1 

Oj  +  2o^02 

2(1  -  V) 

(Oj  +  Oj)^ 

2 

1 

"2 

2(1  -  V) 

(Oi  + 

1 

2vO{ 

2(1  -  v)  Oj  +  Oj  ’ 

1 

2vo2 

2(1  -  v)  Oj  +  Oj  ’ 

2 

1 

2(1  -  V) 

(o,  +  Oj)* 

+  (1  -  2v)- 


ai  +  fl2 


a,  + 


-  (1  -  2v)- 


«i  +  «2 


a,  +  02 


(A-33) 


’1212 


*2323 


1 

f 

2  2 

O,  +  02 

1  -  2v 

2(1  -  V) 

2(Oj  +  Oj)* 

2 

> 

«1 

«2 

2(Oj  +  Oj) 


'3131 


2(Oi  +  Oj) 


It  should  be  emphasized  that  v  in  (A-32)  and  (A-33)  pertains  to  the  Poisson’s  ratio 
of  the  matrhc  material. 
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APPENDEXB 


INCREMENTAL  THEORIES  OF  PLASnCTIY 


B.1  Introduction 

The  purpose  of  this  appendix  is  to  summarize  the  incremental  plasticity 
theories  used  to  calculate  the  increment  of  matrix  plastic  strain,  ACyP™.  Two  theories 
were  utilized:  a  nonlinear  kinematic  hardening  theory  with  backstress  decomposition 
and  an  incremental  isotropic  hardening  theory. 


B.2  Nonlinear  Kinematic  Hardening  Theory 

The  nonlinear  kinematic  hardening  theory  selected  in  this  work  for 
representation  of  matrix  behavior  is  along  the  lines  of  that  employed  by  Chaboche 
and  co-workers  (1978, 1983a,  1983b)  and  further  developed  by  McDowell  (1985)  and 
McDowell  and  Moosbrugger  (1989)  for  nonproportional  cyclic  plasticity.  The 
kinematical  hardening  variable  or  backstress,  Oy,  is  decomposed  into  an  arbitrary 
number  of  components  representative  of  dislocation  interaction  at  corresponding  size 
scales.  This  permits  accurate  description  of  complex  plastici^  phenomena.  Defining 
the  yield  surface  in  the  von  Mises  form 


^ = |W  - «.)(»;'  - «»)  -  (»•« 

where-  =  ajj®  -  (l/3)auj“5ij  and  ajj  is  assumed  to  be  deviatoric,  we  can  write  the 
flow  rule  as 


ej"  =  «»«)”(/  */^=0  ond 

=  0  otherwise. 

By  the  consistency  condition,  if  F=0  during  plastic  flow,  we  may  write 


(B-2) 


/  =  0  =  — dj  ♦  — a„  -  iRk 
So-  '  * 

The  nonlinear  kinematic  hardening  rule  is  of  the  form 


(B-3) 


«<;>  .  C<^b<^nii  -  a'J’Jp  (B-4) 

where  n  =  1,2,...,M  is  the  total  number  of  decomposed  backstresses,  i.e. 


E* 

••1 


i(«) 


(B-5) 


and 


1 


In  equation  (B-2)  and  (B-4), 

I 
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and  and  are  material  constants.  In  this  study,  we  assume  pure  kinematic 
hardening  (R=0).  Using  (B-1),  (B-4),  and  (B-S),  (B-3)  leads  to 


k  =  f; 


(B-7) 


For  the  aluminum  in  this  study,  M=2,  =  12,  =  1400,  =  15,  b^^^ 

=  5.6,  and  R  =  2.0  ksi. 


B  J  Incremental  Isotropic  Hardening  Theory 

Weng  (1988)  tested  his  secant  modulus  theory  on  a  composite  with  an  epo;^^ 
resin  matrix  reinforced  with  silica  particles.  He  gives  the  following  information  to 
represent  the  behavior  of  the  epoxy  matrix: 


o'"  *  o"  +  H(en" 

where  Oy®  =  11  ksi  s  uniaxial  yield  strength;  H  =  4.67  ksi;  n  =  0.26; 


(B-8) 


the  effective  matrix  stress  is 


effective  matrix  plastic  strain  is 


2  ;wi  ;>m 


and  the  deviatoric  matrix  stress  is 


o 


(B-IO) 


(B-ll) 


We  use  isotropic  hardening  plasticity  with  the  yield  surface  (i.e.  no  backstress) 


F  = 


3  m' 

-0,0, 


-  n. 


The  flow  rule  for  determining  plastic  matrix  strain  is 


ej"  =  7(Ojan«)»H  if  F=0  and 
h 

=  0  otherwise. 


where 


(B.12) 


(B.13) 


(B-14) 


The  yield  surface  radius,  Y,  is 


where 


(B-IS) 


» 


» 


||p  -  /||(‘r‘r)^*  =  =  j'"- 

_  # 

This  theory  is  valid  only  for  proportional  loading.  By  the  consistency  condition,  /’=0 

during  plastic  flow;  therefore,  h,  the  plastic  hardening  modulus,  is  given  by 


2  dS* 

■ 


(B-17) 


i 


B-97 


APPENDIX  C 


FLOWCHARTS  FOR  COMPUTER  PROGRAMS 


The  purpose  of  this  appendix  is  to  outline  the  logic  in  the  various  FORTRAN 
programs  employed  to  calculate  the  theoretical  results  presented  in  this  thesis.  This 
appendix  includes  flowcharts  for  the  three  analytical  methods  presented  in  Chapters 
n  and  III.  These  were  the  EIAS  Method  with  Elastic  Constraint  (Program 
BALNIT);  the  EIAS  Method-Tangent  Stif&iess  Formulation  (Program  TANl);  and 
the  EIAS  Method-Tangent  Stiffiiess  Formulation  with  6  parameter  and  x  function 
(Program  CSTAR). 

The  main  program  reads  in  the  stress  rate  and  time  step  from  an  input  file 
called  BALIN.  The  applied  stress  increment  is  calculated  by  two  subroutines  called 
SETUP  and  INTEG  (for  IN  l  EGrationl  The  SCM  (for  Self-Consistent  Method) 
subroutine  is  then  called  to  calculate  the  output  variables.  These  output  variables 
were: 

the  composite  stress; 
the  composite  strain; 
plastic  matrix  strain; 


matrix  stress; 


and  increment  of  plastic  matrix 

strain. 

After  returning  to  the  main  program,  the  stress  and  strain  state  for  that  particular 
time  step  is  written  to  an  output  file  called  BALOUT.  The  program  continues  this 
cycle  until  the  loading  history  is  exhausted. 

The  SCM  subroutine  varies  for  each  of  the  three  programs.  In  the  BALNIT 
program,  the  SCM  subroutine  is  based  on  the  EIAS  Method  with  Elastic  Constraint 
detailed  in  Chapter  II  of  the  main  thesis.  The  flowchart  for  this  SCM  subroutine  is 
shown  in  Figure  C-1.  The  SCM  subroutine  for  the  TANl  program  is  based  on  the 
EIAS  Method  Tangent  Stiffiiess  formulation  presented  in  §  3.1,  and  its  flowchart  is 
shown  in  Figure  C-2.  Finally,  the  SCM  subroutine  for  the  CSTAR  program,  which 
incorporates  the  6  parameter  and  %  function  into  the  EIAS  Method  Tangent  Stiffness 
Formulation,  is  shown  in  the  flowchart  in  Figure  C-3. 


C-99 


From  the  Main  Program:  BALNIT 

Call  SCM  (INIT)  Return  to  the  Main  Program 


Figure  C-1.  Flowchart  for  SCM  subroutine  in  Program  BALNIT. 


From  the  Mein  Program:  CSTAR 

Call  SCM  (INIT)  Return  to  the  Main  Program 


Figure  C-3,  Flowchart  for  SCM  subroutine  in  Program  CSTAR. 


APPENDIX  D 


COMPARISON  OF  THE  TANGENT  STIFFNESS  APPROACH 
WITH  THE  SECANT  STIFFNESS  APPROACH 


According  to  the  Eshelby  approach,  secant  stifbiess  method  (Weng,  1988), 


~  »  (D-1) 

where  =  e^®,  the  comparison  material  strain;  =  ey*,  the  image  strain, 

*  the  perturbation  strain;  and  e^*  is  the  eigenstrain.  Here,  Cp*  is  the  matrix 
secant  modulus. 

Now,  at  the  next  stage  of  deformation,  we  again  apply  the  secant  stiffness 
method,  i.e.. 


Subtracting  state  (D-1)  from  state  (D-2)  gives 

Ci^(Eey  =  iCilEelt  -  tu)  *  -  6ey 

~  ®®h) 


(D-2) 


(P-3) 


C^{Ee||)  ■  -  CjiJj  (D-4) 

*  -  *e-„]  . 

But,  noting  that  the  tangent  stifhiess  approach  is  given  by 

“  **«]  =  ~  ®«)]  » 

where  is  the  tangent  stiffness  of  the  matrix,  it  is  clear  that  (D-4)  may  be  written 
as 


tangent  stifhiess  approach 

Equation  (D-tf)  infers  that  the  tangent  stif&iess  approach,  in  principle,  differs  from 
successive  application  of  the  secant  stiffness  approach  by  only  a  second  order  term 
which  can  be  taken  as  arbitrarily  small  for  small  increments,  provided  the  loading  is 
proportional.  Hence,  the  tangent  stiffness  approach  is  the  logical  analog  of  the 
secant  stiffness  method  for  incremental  matrix  plastici^. 

A  subtle  difference,  however,  is  introduced  if  the  assumption  is  made  in  the 
tangent  stiffness  method  that  the  stress-strain  relationship  of  the  comparison  matrix 
material  is  given  by 


d;  - 


(D.7) 


such  that  =  Qju^  is  assumed.  By  differencing  the  results  of  parallel  solutions 
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from  the  current  state  employing  and  respectively,  in  (D-7),  it  can  be 
shown  that  the  error,  aE,  of  the  tangent  stiffness  approach  with  assumption  (D-7) 
relative  to  successive  application  of  the  secant  stiffness  approach  is  given  by 


a£ 


~  ®</)l 


+ 


(D.8) 


The  notation  0[*]  denotes  that  the  error  is  "on  the  order  of  the  quantity  enclosed 
in  the  brackets.  Qearly,  the  error  is  second  order  provided  that  either 
which  is  highly  improbable,  or  that  ■  Cy^^®,  which  is  perhaps  obvious.  Hence, 
the  assumption  Cyu^°  =  Cy^^  in  the  tangent  stiffness  approach  will  lead  to  first  order 
error  in  representing  the  secant  stiffness  results  unless 


l(ci,  -  C')s{c^-‘a’.  -  C^-'o;)|  0>-») 

is  of  second  order. 

There  are  several  important  regimes  where  (D-9)  is  satisfied,  including  elastic 
behavior  of  both  phases  and  low  volume  fraction  of  reinforcement  or  phase  property 
mismatch  such  that  the  matrix  stress  differs  little  from  the  average  composite  stress. 
Perhaps  the  most  prevalent  and  important  case  is  where  the  matrix  work  hardening 
behavior  is  slight  which  is  common  for  ductile  matrix  materials.  In  this  case,  the 
variation  of  the  magnitude  of  Cyu^  is  small  acro  s  a  wide  range  of  matrix  plastic 
strain. 
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It  must  be  emphasized  that  in  a  fully  incremental  approach,  the  matrix 
behavior  is  permitted  to  be  path  history  dependent,  unlike  in  the  secant  stiffness 
approach.  Precise  assignment  of  is  therefore  an  impossible  task  in  the  general 
case.  This  is  the  basis  for  introducing  the  parameter  B  in  the  tangent  stiffness 
approach,  e.g.  equation  (3-14). 

For  situations  where  the  constraint  of  the  fibers  on  the  matrix  is  high  enough 
in  the  plastic  regime  to  result  in  lay®!  >>  Hffy™!!,  B  is  introduced  to  effectively  force 
the  tangent  stiffness  toward  values  less  than  the  actual  matrix  plastic  stiffness,  even 
in  the  transient  yielding  regime.  This  effect  is  perhaps  most  relevant  in  off-axis  and 
shear  loadings  since  matrix  plasticity  has  little  influence  on  behavior  for  loading 
nearly  aligned  with  the  fiber  direction. 

In  summary,  it  appears  that  the  tangent  stiffiiess  model  in  its  general  form  is 
a  quite  simple,  versatile  and  accurate  approach  for  incremental  elastic-plastic 
composite  deformation.  It  requires  no  iteration  at  each  loading  increment  and 
requires  only  constituent  properties  for  user  input. 
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APPENDIX  E 


CONSTITUENT  PROPERTIES 


I 

E.1  Introduction 

We  investigated  the  inelastic  behavior  of  two  composite  systems.  The  boron- 
'  aluminum  system  is  a  unidirectional,  continuous  fiber-reinforced  composite  for  which 

the  fibers  are  boron  and  the  matrix  is  aluminum.  The  silica-epoxy  system  is  a 
spherical  particle-reinforced  composite  where  the  spherical  particles  are  made  of 
silica  and  the  matrix  is  an  epoxy  resin.  Listed  below  are  the  elastic  properties  of 
these  constituent  materials  which  were  used  in  our  investigation. 

I 

E.2  Constituent  Properties  of  the  Boron-Aluminum  System 

The  boron  fibers  are  considered  elastic  while  the  aluminum  matrix  is  elastic- 

» 

plastic  in  the  B/Al  system.  Both  the  aluminum  matrix  and  boron  fibers  were  assumed 
isotropic.  Their  respective  elastic  properties  are  shown  below. 

» 


) 


r 


■ 


Ea 

Ga 

Va 

Er 

Gj 

(lO^psi) 

(lO^psi) 

(lO^psi) 

(lO^psi) 

[GPa] 

[GPa] 

[GPa] 

[GPa] 

6061 A1  10.5 

3.95 

0.33 

10.5 

3.95 

[72.5] 

[27.2] 

[72.5] 

[27.2] 

Boron  58.0 

24.2 

0.20 

58.0 

24.2 

[400.0] 

[166.8] 

[400.0] 

[166.8] 

where  the  subscript  A  corresponds  to  the  33  (fiber)  or  the  "aligned"  direction  and 
the  subscript  T  corresponds  to  the  12  or  the  "transverse"  plane.  E,  G,  and  v  are  the 
Young’s  modulus,  shear  modulus,  and  elastic  Poisson’s  ratios,  respectively.  These 
properties  were  obtained  from  Dvorak  and  Bahei-El-Din  (1987).  The  material 
properties  used  to  characterize  the  flow  of  the  matrix  are  listed  in  Appendix  B. 

E3  Constituent  Properties  of  the  SUica<Epoi^  System 

In  the  silica-epojqr  system,  the  silica  is  elastic  and  the  epo^g^  is  elastic-plastic. 
Both  phases  are  assumed  to  be  isotropic.  From  Weng  (1988),  the  elastic  properties 


are  as  follows: 


E 

(lo^psi) 

[GPa] 

V 

Silica 

10.6 

0.18 

[73.1] 

Epoxy 

0.458 

0.35 

[3.16] 

The  material  properties  used  to  characterize  the  flow  of  the  matrix  are  listed  in 
Appendix  B. 
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APPENDIX  F 


ANISOTROPIC  CONSTRAINT  HARDENING 
AND  THE  X  FUNCTION 


F.l  Introduction 

In  this  thesis,  we  conducted  a  theoretical  analysis  of  the  inelastic  behavior  of 
metal-matrix  composites  reinforced  with  aligned,  continuous  fibers  and  attempted  to 
model  the  experimental  response  of  the  Boron-Aluminum  system.  Due  to  the 
anisotropy  associated  with  the  aligned  nature  of  the  reinforcement  and  the  high 
transverse  shear  stiffness  of  the  boron  fiber,  constraint  hardening  is  highly  anisotropic 
in  the  boron-aluminum  composite  system. 

In  our  tangent  stiffness  model  (§  3.1),  we  observed  a  hardening  response  in 
the  plastic  regime  for  off-axis  loading  because  we  did  not  take  this  constraint 
hardening  phenomenon  into  consideration.  Further  investigation  revealed  that  relief 
of  constraint  hardening  in  this  composite  system  was  due  to  rotation  of  the  matrix 
plastic  strain  rate  vector  away  from  the  fiber  direction  towards  the  transverse 
direction  (which  is  the  direction  of  least  constraint).  We  incorporated  this  mechanism 
for  constraint  relief  into  our  model  via  the  %  function.  With  this  %  function  we 
computed  an  effective  fiber  stiffness  which  ensured  a  continuous  transition  to  matrix- 
dominated  behavior  as  is  observed  in  the  plastic  regime  for  the  boron-aluminum 


composite  system. 


The  purpose  of  this  appendix  is  to  explain  this  phenomenon  of  constraint 
hardening  as  it  relates  to  composites  reinforced  by  aligned,  continuous  fibers  and 
provide  a  justification  for  the  form  of  our  proposed  %  function.  Furthermore,  we  will 
discuss  how  it  might  apply  to  different  types  of  composite  systems. 

F.2  Anisotropic  Constraint  Hardening 

To  understand  constraint  hardening,  let  us  consider  an  element  of  matrix 
material  which  is  bounded  on  two  sides  by  fibers  (see  Figure  F-1).  For  an  applied 
stress  in  the  X3-direction,  the  matrix  material  cannot  contract  laterally  in  the  Xj- 
direction  due  to  the  high  transverse  shear  stiffness  of  the  boron  fibers.  This  lateral 
constraint  results  in  a  stiffer  response  for  loads  applied  in  the  fiber  (X3)  direction. 
However,  when  a  stress  is  applied  in  the  transverse  (Xj)  direction,  the  matrix  deforms 
"naturally"  or  in  a  matrix-dominated  mode  since  there  is  no  constraining  effect  on  the 
top  and  bottom  surfaces  of  the  element.  We  can  see  that  the  level  of  constraint 
varies  considerably  betwee  *he  fiber  and  transverse  directions.  Therefore,  we  say 
that  the  constraint  hardening  is  anisotropic.  It  is  important  to  realize  that  the  fiber 
direction  is  the  direction  of  greatest  constraint  while  the  transverse  direction  is  the 
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direction  of  least  constraint. 


F3  Relief  Mecbanism  for  Anisotropic  Constraint 

It  has  been  experimentally  observed  for  off-axis  loadings  that  boron-aluminum 
exhibits  a  transition  to  the  so-called  "matrix-dominated  mode"  of  deformation  in  the 
plastic  regime.  Therefore,  it  is  evident  that  there  is  a  relief  of  constraint  hardening 
effects  associated  with  the  presence  of  the  stiff  boron  fibers.  We  must  identify  the 
mechanism  responsible  for  this  relief  of  constraint  in  order  to  successfully  model  the 
experimental  behavior. 

Dvorak  and  Bahei-El-Din  (1988)  observed  that  even  for  proportionally  applied 
off-axis  loadings  there  is  nonproportionality  of  the  plastic  strains  in  the  matrix.  For 
an  off-axis  loading,  the  plastic  strain  rate  vector,  Cy**®,  is  initially  collinear  with  the 
matrix  stress  rate  (see  Figure  F-2).  However,  as  plastic  flow  develops,  ey**®  rotates 
away  from  the  fiber  direction  towards  the  direction  of  least  constraint.  This  rotation 
of  CyP®  accounts  for  the  observed  nonproportionality  of  the  plastic  matrix  strains.  At 
this  point,  there  is  sufficient  rationale  to  conclude  that  the  rotation  of  Sy^  is 
responsible  for  the  observed  relief  of  constraint  hardening.  To  confirm  this  rationale, 
we  analyzed  the  rotation  of  the  matrix  plastic  strain  rate  vector,  the  results  of  which 
were  presented  in  §  4.5  of  the  main  thesis.  We  will  restate  those  results  here. 

We  utilized  the  results  from  our  tangent  stiffness  model  to  analyze  the  rotation 
of  eyP®.  It  is  important  to  remember  that  the  tangent  stiffness  model  accurately 
approximated  the  response  of  the  composite  in  the  initial  yielding  regime.  Therefore, 
it  is  accurate  up  to  approximately  0.1  or  0.2  %  strain.  Hence,  we  are  justified  in 
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using  the  results  of  this  model  in  our  analysis  of  eyP“. 

Referring  to  Figure  8,  we  see  that  in  the  0°  orientation  the  plastic  strains 
remain  proportional  as  indicated  by  the  linearity  of  the  plot  of  633*"”  versus  £22^' 
However,  in  the  off-axis  orientations,  there  is  curvature  of  the  plot  indicating 
increasing  dominance  of  622'’“  and  rotation  of  6^**“  away  from  the  fiber  direction 
towards  the  direction  of  least  constraint,  the  transverse  direction.  This  analysis 
confirms  our  hypothesis  that  rotation  of  ey*””  is  the  mechanism  responsible  for  the 
relief  of  constraint  hardening  and  the  transition  to  a  matrix-dominated  response  in 
the  boron-aluminum  composite  system. 


F.4  Analyzing  Constraint  Hardening  Utilizing  Wang’s  Representation  Theorem 
At  this  point,  we  have  established  that  there  is  decreasing  constraint  hardening 
in  the  boron-aluminum  composite  associated  with  the  development  of  plastic  flow  in 
the  matrix  due  to  rotation  of  eyP“.  K  we  define  i|ry  as 


then  i^y  is  clearly  the  rate  of  change  of  the  hydrostatic  constraint  stress,  with 
plastic  matrix  strain.  In  general,  we  may  represent  tjry  as  an  isotropic  tensor  function 
of  the  form 


(F-2) 
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The  unit  vector  63  represents  constraint  hardening  due  to  the  anisotropy  in  the 
composite  due  to  fiber  reinforcement  while  represents  constraint  hardening  in 
the  matrix  due  to  the  impedance  of  dislocation  motion  by  the  presence  of  the 
reinforcement.  The  semicolon  indicates  implicit  dependence  of  ijr^  on  the  functions 
B  and  x  which  are  our  proposed  methods  of  representing  constraint  hardening  effects 
in  the  composite.  Using  Wang’s  representation  theorem,  i|rjj  may  be  expressed  in 
terms  of  the  irreducible  functionality  basis  of  the  pair  (ej.cry'”)  (Wang,  1970)  as 
follows: 

(F-3) 

Here,  9^  is  expressed  as  the  sum  of  four  isotropic  tensor  functions  whose  arguments 
are  as  shown  above.  The  first  term  in  the  expression,  cpi,  is  included  based  on  the 
formality  of  the  representation  but  will  drop  out  since  the  first  invariant  of  is 
equal  to  zero.  The  second  term,  (pj,  represents  constraint  hardening  in  fiber- 
reinforced  composites.  The  third  and  fourth  terms,  93  and  9^,  represent  constraint 
hardening  in  particle-reinforced  composites.  Since  the  second  term  and  the  third  and 
fourth  terms  are  representative  of  two  completely  different  classes  of  composite 
materials,  cross  terms  in  the  representation  have  been  omitted.  Also,  joint  invariants 
of  Ou“  have  been  omitted  as  second  order. 

We  will  now  take  a  closer  look  at  constraint  hardening  in  fiber-reinforced 
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composites  using  Wang’s  representation  theorem.  For  fiber-reinforced  composites, 
we  may  in  general  express  tlry  using  the  second  term  in  (F-3)  as 

From  §  3,2  of  the  main  thesis,  we  recall  that  constraint  hardening  due  to  is 
isotropic  and  that  constraint  hardening  associated  with  S3  is  anisotropic.  We  may 
think  of  6  as  the  function  utilized  to  represent  isotropic  constraint  hardening  while 
X  is  the  function  utilized  to  represent  anisotropic  constraint  hardening.  If  we  refer 
to  Figure  F-3,  we  can  see  that  the  anisotropic  constraint  hardening  due  to  fi  and  % 
is  dominant  over  isotropic  constraint  hardening  for  fiber-reinforced  composites.  This 
graphical  illustration  in  combination  with  our  analysis  in  §  3JI  leads  us  to  conclude 
that  isotropic  constraint  hardening  is  of  second  order  importance  in  the  fiber- 
reinforced  composite. 

F,5  Development  of  the  %  Function 

We  have  shown  that  relief  of  constraint  in  the  fiber-reinforced  composite  is 
associated  with  the  rotation  of  ejj*®*  of  away  firom  the  fiber  direction  during  off-axis 
loadings.  Furthermore,  we  have  shown  that  constraint  hardening  is  an  anisotropic 
phenomenon  in  this  class  of  composites. 

'  From  §  3,2,  the  hydrostatic  constraint  stress  rate  is 

-  siCp-b  -  ,  (F-S) 
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where  n33  is  the  component  of  n^j,  the  unit  normal  vector  in  the  direction  of  eij*™,  in 
the  fiber  direction,  63  and  %  is  an  isotropic  scaler  function  of  63  and  As  plastic 
flow  in  the  matrix  develops  and  6^^“  rotates  away  from  63,  n33  decreases  in 
magnitude.  We  also  know  from  experimental  evidence  that  the  boron-aluminum 
system  demonstrates  a  rapid  transition  to  matrix-dominated  behavior  (see 
experimental  data  plotted  in  Figure  10,  for  example). 

Utilizing  these  facts,  we  are  able  to  model  constraint  relief  by  effectively 
reducing  the  stiffness  of  the  fiber  utilizing  the  %  function  as  follows; 

where  %  takes  on  the  form 


The  idea  here  is  to  facilitate  the  transition  to  matrix-dominated  behavior  (the  tangent 
stiffness,  CyJ)  by  effectively  reducing  the  fiber  stiffiiess  from  Cgu*  in  the  elastic 
regime  to  in  the  inelastic  regime. 

Some  explanation  of  the  terms  in  (F-7)  is  required.  The  asymptotic  plastic 
hardening  modulus,  h^,  is  indicative  of  the  asymptotic  plastic  behavior  of  the  matrix 
material.  Therefore,  the  ratio  h^  is  representative  of  the  plastic  strain  state  of  the 
matrix.  In  the  elastic  regime,  where  h  -» «,  h^/h  -*  0  and  x  1.  However,  as  plastic 
flow  develops,  ho/h  1  and  contributes  to  the  overall  value  within  the  brackets. 
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The  component  of  Hjj  in  the  63  direction,  n33,  is  indicative  of  the  rotation  of 
ejjP“.  As  plastic  flow  develops  and  constraint  is  relieved,  njj  -►  0.  As  indicated  by  the 
experimental  data,  matrix-dominated  behavior  is  realized  earlier  on  for  more  off-axis 
orientations  (  e.g.  45®)  than  for  off-axis  orientations  closer  to  the  fiber  direction  (e.g. 
10®).  This  is  because  n33  is  initially  greater  for  the  10°  orientation  say  than  the  45® 
orientation.  Therefore,  it  takes  longer  for  matrix-dominated  behavior  to  develop  in 
the  10®  orientation. 

As  stated  above,  the  boron-aluminum  system  exhibits  a  rapid  transition  to 
matrix-dominated  behavior.  This  may  not  be  the  case  for  other  composite  systems. 
The  exponent  N  allows  us  to  dictate  the  rate  of  transition  to  matrix-dominated 
behavior.  In  general,  the  higher  the  transverse  shear  stiffness  of  the  fiber,  the  faster 
the  transition.  For  the  boron-aluminum  system,  we  fitted  to  the  15°  orientation  using 
N=4. 

The  aspect  ratio  a/c  is  included  to  take  into  account  the  nature  of  the 
reinforcement.  For  aligned,  continuous  fibers,  a/c  =  0.  For  spherical,  particle 
reinforcement,  a/c  =  1. 

Taking  the  above  factors  into  consideration,  we  can  see  that  the  %  function 
generally  behaves  as  shown  in  Figure  F-4.  When  there  is  no  plastic  flow,  x  =  1  and 
Ajki*  ~  Qjki*  all  off-axis  orientations.  As  plastic  flow  develops,  x  exponentially 
decays  to  zero  at  a  rate  which  is  dependent  on  the  loading  orientadon  and  the 
transverse  shear  stiffness  of  the  fiber.  Figure  F-4  is  intended  to  give  a  general 
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representation  of  this  behavior  and  is  not  indicative  of  a  particular  system. 

F.6  Conclusion 

In  summary,  the  %  function  is  a  mathematical  means  of  modeling  the  relief  of 
anisotropic  constraint  hardening  in  fiber-reinforced  composites.  It  is  theoretically 
based  on  the  fact  that  constraint  is  relieved  by  rotation  of  away  from  63  during 
the  development  of  plastic  flow  of  the  matrix  and  that  matrix-dominated  behavior  is 
exhibited  in  the  plastic  regime. 

Its  development  was  based  solely  on  our  attempts  to  model  the  inelastic 
behavior  of  boron-aluminum.  Applications  to  other  systems  such  as  graphite- 
aluminum  have  as  yet  to  be  investigated.  Isotropic  constraint  hardening  in  the  boron- 
aluminum  composite  was  neglected  as  second  order  because  of  the  dominance  of 
anisotropic  constraint  due  to  the  high  shear  transverse  sti&ess  and  anisotropy  of  the 
reinforcement.  Therefore,  the  6  parameter  and  the  concept  of  the  comparison 
matrix  material  were  not  incorporated  into  our  attempts  to  model  boron-aluminum. 
We  recognize  that  for  other  systems,  this  may  indeed  not  be  the  case.  A  combination 
of  B  and  X  nt^y  be  required  to  successfully  model  the  inelastic  behavior  of  systems 
with  more  compliant  reinforcement. 
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Figure  F-2.  Rotation  of  the  matrix  plastic  strain  rate  vector  away  from  the  fiber 
direction  during  the  development  of  plastic  flow  in  the  matrix  material. 
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Figure  F-3.  Anisotropic  constraint  hardening  in  fiber-reinforced  composites. 
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Figure  F-4.  Rate  of  decay  of  %  for  various  off-axis  orientations  as  a  function  of 
plastic  matrix  strain. 


